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Abstract 


A synthetic aperture acoustic array is used to as a spatial filter in a reverberant test 
facility. The array is used to examine the c_ snge in the received signal when a fluid 
loaded steel plate is replaced by the free surface. Twenty millisecond active pulses at 2 
kHz were scattered by the target plate and recorded by the synthetic aperture array for 
post processing. T 1e method of reverberant field subtraction was utilized and found 
unsatisfactory for use with signals of 20 msec duration in a reverberant tank. Array 


processing using a rectangular amplitude function was marginally successful. A 


Chebyshev amplitude weighted array was successfully used to demonstrate the ability of 


a linear array to discriminate between the steel plate and free surface scattering. 


A computer aided interactive design tool was written to aid in the analysis of array 
architectures by calculating hydrophone amplitude weights, beamwidth, directivity index 
and graphically presenting the normalized directional beampattem. The program allows 


comparative analysis of two array designs simultaneously. 
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Nomenclature 


a(t) Amplitude of the transmitted periodic signal 

c Sound speed in feet/second 

d Distance between adjacent hydrophones 

DI Directivity index 

F(O) Array factor 

G(9) Normalized directional gain pattem 

iF Amplitude weight of j" hydrophone 

1, Modified Bessel function of order 0 

k Wavenumber 

m,,m, Mean of temporal signals x(t), y(t) 

M Number of elements constituting an array 

n Taylor parameter to determine a.ceptable sidelobe levels 
p Normalized array length (-1.0Sp<1.0) 

P Pressure 

R,, Cross correlation between x(t) and y(t) 

Rap Desired sidelobe level in dB 

Sp(t) Signal at the hydrophone 

S7(t) Signal at the projector 

t, Pulse duration 

Zo Chebyshev polynomial 

a Pedestal height for cosine amplitude taper function 
, Variable used in calculation of amplitude tapered array factor 


Wavelength in feet 











Carrier frequency in radians/second 

Phase modulation 

Correlation coefficient or normalized covariance 

Variance of temporal signal x(t) 

Time delay between signal arriving at adjacent hydrophones 
Angular direction from array axis 


Main response axis, “look direction" 




















1 INTRODUCTION 


1.1 Purpose of the Research 

The scattering of acoustic energy from fluid loaded structures has been a topic of 
interest for many years. Recently attention has shifted to considering active control 
methods to reduce/eliminate the energy reflected from these structures. An active control 
system would be capable of detecting the incoming wavefront, determine its physical 
properties, and through an adaptive control system cause the structure to respond ina 
manner which will reduce the reflection of the incident wavefront. 

As these studies mature, the need for inwater experiments becomes necessary to 
validate the theory and test the control system algorithms. While experiments in open 
water would be most like the anticipated operating environment they would also be 
prohibitively expensive. A likely alternative is to conduct the inwater tests in a tank 
facility. This approach has a problem in that the physical size of the tank lends itself to 
developing a large number of reflections of the original signal off the sides, botton:, and 
free surface tending to corrupt that part of the signal which undergoes specular reflection 
from the fluid loaded structure. The researcher requires quantifying the change in 
scattered energy from the structure due to the control system in order to rate the 
performance of a chosen method. 

A simple solution to the problem is to arrange the experiment in the tank such that 
the scattered signal reaches the hydrophone before the second multipath signal. In this 
way time clipping the signal appropriately will show the seflected signal alone. Figure 1 
shows a typical setup for this type of experiment. 

The significant drawback to this method is that it artificially requires the active 


interrogation of the structure to be a short duration pulse in order to assure that the 


desired scattered signal is discernable in the hydrophone output. As the associated signal 
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Figure 1.1 Timeclipping to eliminate multipath interference 


time history implies, separation of the desired signal from direct path and multipath 
signals may be difficult. The short active pulse reduces the time available to the control 
system for detection, determination of the signal characteristics, and initiation of an 
appropriate control sequence. Alternatively, and the purpose of this report, a series of 
hydrophones can be used to create an array, a spatial filter used to discriminate signals 
from directions other than the desired look vector. Figure 1.2 shows the previous 
experiment with a 9 element array used to spatially filter out the signals from non-look 
directions. The main lobe of the array is directed at the target of interest while the 
projector is placed close to a null in the beam pattern minimizing detection of the direct 
signal. 

The use of a tank for the experiment creates an undesirable environment for an 


array. Multipath signals are numerous and correlated. Phase cancellation of the signal 


1] 

















ACTIVE 
CONTROL 
SYSTEM 


SHAKER 


PROJECTOR 


@ 


LINEAR ARRAY 
OF NINE 
HY DROPHONE S 










" 

C 

= 

= YC 

4 SUMMED QUTPUT 
Oo SIGNAL 

Zz 


Figure 1.2 Nine element array used as a spatial filter 


from the desired direction is a strong possibility. An advantage of array processing is its 
ability to filter signals from undesired directions allowing longer active interrogation of 


the structure which in turn gives the control system more time to respond. 


1.2 Research Goals 


The purpose of this research is to evaluate array processing for use in supporting 
continued experiments in the active control of scattering. The objective is to determine a 
beamforming scheme supporting a longer active interrogation of the structure. This in 
turn gives more response time to the detection and active control system. 

As currently envisioned the initial experiment in active scattering will be run using 


a time clipped short duration pulse of approximately 2 msec to eliminate any multipath 
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signal detection. The goal of this research is to determine an array beamforming method 


which will give an order of magnitude increase in the active pulse length. The stated 
objective is to allow active interrogation of up to 20 msec. This increase should greatly 
enhance the quantitative assessment of the active control system performance. 

A secondary deliverable product as a result of this research is an interactive 
computer aided design tool which can be used to assist in visualizing the spatial gain due 
to array processing. The program models the array gain of two separate designs such that 
a user could easily discriminate between the tradeoffs associated with each. It is intended 
that this program be used as a first estimation of the characteristic design parameters of a 


linear discrete element array. 


1.3 Organization and Content 


The report is divided into four major sections; theoretical analysis experimental 
setup, computer support, and analysis of results. The first, theoretical analysis, discusses 
the mathematical development of the beamforming algorithms used as well as an estimate 
of their applicability to the specifics of beamforming in a reverberant tank. The 
numerical methods employed in creating an interactive program used to design a 
multielement line array are presented. The major subheadings include beamforming, 
phase and amplitude taper, and output power calculations. The second section deals with 
the experimental setup in the reverberant acoustic facility as well as methods of data 
collection, reduction, and processing. The third section deals with an interactive array 
design tool written to assist in determining the capabilities of various array architectures. 
The final section deals with an analysis of the experimental results as well as 


recommendations for future study. 








2 Theoretical Analysis 


2.1 Signal Characterization 


The following discussion and throughout this paper the nomenclature of Elliott 
[15] is used. The transmitted signal, a finite duration pulse in this study, is denoted by 


S,(t) at the transmitting element (projector). 


a(t) cos(wt + 0(1)) O<rs< “| 


he 2.1 
st) 0 everywhere else on 


Where a(t) is the amplitude, @ the carrier frequency, @(t) the phase modulation, and 
t, the pulse duration. !t is also possible to examine the signal at a receiving element 


(hydrophone). For a hydrophone located a distance r from the projector the received 


sp(t) ataft 1 (2.2) 


Noting that c is the wave propagation speed in water it is obvious that the received signal 


signal would be of the form: 


is proportional to the time delayed projected signal and that this time delay is simply that 
due to the wave progressing from the projector to the hydrophone. It is common practice 
to denote the signals represented by equations 2.1 and 2.2 as equivalent complex signals 
where the actual signal is the real part of the complex representation. 

S p(t) = VP et ™ *#* (2.3) 

= cos(cor + (t) + 6) + sin(wr + (1) + 8) (2.4) 
Where s,(t) is equal to the real part of Equations 2.3 and 2.4. 
It is often convenient to examine the information carrying part of the signal by 

eliminating the carrier. The baseband reduced complex envelope, Pillai [13], or complex 


baseband demodulated version Elliott [15], of a signal eliminates the carrier portion of 
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the signal leaving a real inphase part and an imaginary quadrature part. Creation of the 
baseband reduced complex envelope can be accomplished by using a Finite Impulse 
Response (FIR) linear phase lowpass filter [15]. 


Sa(t) = VP ef*® (2.5) 


=i, (t) +q,,(¢) (2.6) 


2.2 Statistical Considerations 


Acoustic signal processing often relies on a statistical representation of the signals 
encountered. Relying on the definitions and nomenclature of Newland [5], for two time 
varying signals x(t) and y(t) with means m, and m, and variances 6, and o, respectively, 


the correlation coefficient or normalized covariance is: 


_E{(x—m,)(y —m,)} 
~ 0,0, 


xy 


The value of p,, lies between 0 and 1. A value of 1 indicates coherent signals. A value of 
0 indicates uncorrelated signals. In fact, most signals lie somewhere between these 
extremes. 

Defining E[f(t)] as the expected value of the function f(t), the cross correlation 


between two functions of time is defined as: 


R(t) =E[x(y(t +9) 


Both cross correlation and correlation coefficient will be referred to throughout this 


report. 
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2.3 Conventional Beamforming Analysis 


2.3.1 Enhancements Due to Array Processing 





Arrays have become common in both radar and sonar applications. The main 
function of array processing is to create a larger Signal to Noise Ratio (SNR) by using 
judiciously selected time or phase delays which cause coherent addition of the desired 
signal and incoherent addition of noise and interfering signals. Figure 2.1 shows a 


hypothetical source s(t) and a four element linear array of length L and element spacing 


d’. 











\ <2 
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Figure 2.1 Linear array coherent addition 


If the source is sufficiently distant from the array we can assume that the energy 
arrives in the form of an advancing plane wave from a direction © with respect to the 
array axis. It is seen that the wavefront passes array element 1 at time t, and element 4 at 
time t,. Noting again the assumption of plane wave progression the wave travels a 
distance of L cos © at speed c between element 1 and element 4. Simply stated the signal 


at element 4 at time t is the same signal that passed element | some time (7) earlier. 


Letting n refer to the element number and c be the wave celerity: 








LcosO 
t=—— 
Cc 


_n(d')cos© (25) 
c 

If the signal received at hydrophone 1 sp,(t) is known, as well as the source frequency, we 
can express the expected signal due to the source at each successive hydrophone (n) as 


Sp,(t) = Sp,(t-T). 


Spn(t) =Sp)(t —T) (2.6) 

=s.(he" (2.7) 
=Se(ne ‘ (2.8) 
=Sp(ne  * (2.9) 


By phase shifting, hydrophone outputs the signal from a desired “look direction” can be 
added coherently along the array. The actual output of any hydrophone in the array is a 
sum of the signals from all sources as well as other noise. If the noises are considered 
uncorrelated and the other sources are sufficiently separated trom the desired look 
direction in angular space then Pillai [15] and Steinberg [16] demonstrate that for a single 
source case with uncorrelated identical element noise, the array SNR is a factor of M 


greater than singie hydrophone element SNR. M is the number of elements in the array. 


2.3.2 Distributed Array Normalized Directional Pattern 

Figure 2.2 shows a schematic representation of an array beamformer. Each 
hydrophone element is connected to a phase shifter. The phase shifted signal is then 
sumuned over the array to create the output signal z(t). The phase shifters allow the array 


to be spatially selective as explained in the previous section. Shifting the look direction 
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changes the weighting factor for each element and effectively steers the array. 
Continuously changing the weighting factors will cause the array to scan through angular 


space. 
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Figure 2.2 Schematic of a Beamformer 


The directional pattern can be found by analyzing the array factor defined in 
equation 2.10 which is simply the summation of the element time shifts over the M 
elements constituting the array. 


F@)=S 6" (2.10) 


j=l 
Calculation of the directional pattern can be accomplished by examining the summation 
of the finite series. Using equation 2.10 with an interelement spacing of d for the 3 


element array: 


3 cos 
F(@)= Se" 
tal 
= 2 E, cof a cos) 
=1+ cos 2h is36 + cos 2n2~ cos® 
r r 
18 





Recognizing the periodic nature of cosine functions the second two terms can be 


combined. 


F(O)=1 +204 a cos6| 


The array directional pattern is the result of a finite series the number of terms depending 


on the number of array elements. Dyer [17] generalizes for an odd number of elements: 








2 cos aed cos @] sin od cos @] 





. rd cos® 
m a 


And for an even number of elements: 


cos = d cos 0) sin( = d cos 6] 
, ad 


F(@)=2 
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Multiple sources, Pillai [13], Widrow [14], Steinberg [16], derive the simplified 
normalized directional gain pattem for an array of distributed periodic elements. The 
directional or radiation pattern is the same regardless of even of odd number of array 
elements and is of the form in equation 2.12. To visualize the beam pattern the value of 


© is varied over the spatial angle of interest (ie 0 to 180 degrees). 
a 
G(O) = ——_—__ (2.12) 
M sin—— 


Simplifying Equation 2.12 by substituting k=2n/A leaves: 
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CON ines (2.13) 


Steering the beam to a desired look direction is accomplished by adding a constant phase 
factor to each of the element outputs. This has the effect of shifting the main beam and is 
incorporated in equation 2.11 by substituting cos(©) -cos(©,) for cos(©) where ©, is the 


desired steer angle. 
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Figure 2.3 Six element array steered to broadside 


Equations 2.12 through 2.14 are used to make first approximations as to the ability 
of an array to meet certain design objectives. Assuming that the physical properties of 


the signal and the medium (ie frequency, sound speed etc...) are established for the 
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designer only the features of the array are capable of modification. Limiting the 
discussion to omnidirectional elements in a linear pattern allows the array designer to 
vary the number of elements, interelement spacing, and steer direction to create an 
acceptable array. 

Figures 2.3 through 2.7 are plots of normalized directional patterns using equations 
2.13 and 2.14. In all cases the array elements are considered omnidirectional 
hydrophones, sound speed is 4800 feet/second, the signal is a 2 kHz sin wave (A = 2.4 
feet), and the arrays are linear. 

Figure 2.3 shows a six element array steered to broadside, perpendicular to the array 
axis. The interelement spacing is one half wavelength. The directional pattern is plotted 
as a log function while the direction along the x-axis is linear. In the coordinate system 
established in Figure 2.1 the steer direction angle is from the array axis. Other sources, 
Steinberg [16], Dyer [17], like to define the look angle from broadside, which simply 
replaces the cosine function in equations 2.12 through 2.14 with a sine function. 

in Figure 2.4 the number of elements are increased from six to twelve keeping 
interelement distance at one half wavelength. Both the 6 element and twelve element 
arrays are superimposed to show the effect. The side lobes are reduced and the main lobe 
is narrower while the array is twice the size. It is interesting to note the significant nulls 
between the sidelobes and the major null at endfire. 

Figure 2.5 demonstrates an important point, spatial aliasing of the signal. If the 
interelement distance is increased to exactly one wavelength, it is easy to show that the 
signal will nov: add coherently at the endfire angles where it was previously a null. This 
results in a spatial gain at end fire equal to that of the mainlobe, definitely not a desired 
attribute. For this reason interelement spacing should not normally exceed one half 
wavelength. Keep in mind that for any chosen spacing there is a frequency above which 


aliasing will occur. 
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Figure 2.4 Twelve element array steered to broadside 


This effect creates grating lobes which significantly reduce the effectiveness of the 
array to act as a spatial filter. In order to eliminate grating lobes the interelement distance 
should not exceed A/2 for the highest frequency of interest. An alternate visualization is 
that grating lobes will always exist but judicious selection of interelement separation will 
force them to occur outside the visible region (0°SO<360"). 

Decreasing the interelement spacing below one half wavelength will widen the 
beams and increase the sidelobe level. Another simple calculation will verify that a 6 
element array with A/2 spacing will have the same beam pattem as a 12 element array 
with a spacing of A/4. The difference though is that at any higher frequency of signal the 
6 element array will begin to exhibit spatial aliasing while the 12 element array will 


actually improve in perfonnance. 
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Figure 2.5 The effects of interelement spacing on aliasing. 





Figure 2.6 uses equation 2.13 to show how an additional phase shift factor can be 
used to steer the array from broadside to some other desired look direction. Here the 
array is steered to an angle of 45 degrees forward of broadside. This figure clearly shows 
a second main response axis beam at 315°. Linear arrays form beams which are in fact 
cones symmetric about the array axis. The beampattern plots are two dimensional 
sections of the three dimensional beams showing both sections of the cone. Correct 
visualization of the area in physical space that the main beam covers requires rotating the 
beampattern plots through 180° in angular space. Figure 2.7 is a representation of the 3 
dimensional beams formed by a linear array. Beams are shown at endfire, 30°, 60°, and 


broadside tot the array axis. 
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Figure 2.6 Twelve element array steered 45 degrees 











Figure 2.7 Linear array beams in 3-D. 
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2.3.3 Effects of Array Amplitude Taper 


Equation 2.10 which defined the array factor F(©) was simplified in that the 
individual elements were uniformly spaced and evenly weighted with an amplitude taper 
of unity across the array aperture. A more stringent definition is equation 2.15 where it is 
possible to individually weight the amplitude of each element's output. Each element 
now can be phase tapered as well as amplitude tapered to create the most desirable 
directional gain pattern. 


M an cosO 
F(@)= 5 jek (2.15) 


pat 

In an attempt to quantify the merit of an array the designer relies on three factors. 
Main lobe beamwidth is calculated at the half power point or alternately at the point 
where the normalized directional gain, G[8], is reduced by 3 dB. Sidelobe level, the 
second factor, is the value of G[{6] at the peak of the first sidelobes. Directivity factor (d) 
is the expected (average) value of (F[@])’ over angular space and is a dimensionless 
factor. Directivity factor (d) has a value less than 1.0. The lower the value of d the more 
efficient the array can be considered. 


freer 


a at 





The directivity index, DI, is the ratio, in decibels, between the power received by an 
omnidirectional hydrophone in an isotropic field to the power received by the array beam. 
The larger the value of DI the more effective the array in filtering signals from 


non-desired directions. 
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I 
Directivity Index = 10a 1 


47] 
fr@y 


Amplitude tapers have several effects associated with their use. 


= 10log 


* Increased taper can decrease the associated side lobes, perhaps significantly. 
This is the single most important positive attribute of amplitude taper. 

« Amplitude taper decreases the directivity index a small amount. 

¢ Main lobe beamwidth is increased as is the angular distance between zeros. 
The design program ARRAY in Appendix A allows the user to select several typical 
amplitude taper functions. Altematively the user can enter specific values of amplitude 
taper for each element and obtain a graphical representation of the expected beam pattern. 

Incorporation of variable amplitude weighting in equation 2.15, eliminates the 
ability to derive equations 2.13 and 2.14 by successive simplification of the finite series. 
As before, the array factor F(©) varies depending on an even or odd number of array 
elements, but unlike the uniformly tapered case, the final result is unique for each case. 
A generally applicable equation, such as 2.14, is no longer available. Figure 2.8 shows 
two linear arrays of equally spaced elements steered to an arbitrary direction O. The 
reference position for determining hydrophone location is the center of the array. Array 
A has an odd number of elements and Array B has an even number of elements. 

Balanis [18] solves equation 2.15 for the two cases portrayed in figure 2.7 and 
arrives at equations 2.16 and 2.17 for the array of even and odd elements respectively. A 
change in nomenclature occurs here where the array has either 2M or 2M+1 elements 


where M is an integer value. 
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Figure 2.8 Amplitude tapered arrays of even and odd elements 





F(O).y4 = 5 iF cos ar ia cos(0)| (2.16) 
jel 
M+1 J. 
F(©).y41= x 6 COU - kd cos(@)) (2.17) 


It should be noted that «, is defined as: 
2 j=] 
of Jal 
PNY J#1 
There are a multitude of taper functions applicable to signal processing. The design 
tool ARRAY written for this research allows the user to interactively select from a menu 


of six choices. The numerical methods used in incorporating these techniques are 


explained in detail in Chapter 4. The rectangular, cosine and raised cosine apertures are 
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based on the presentation by Steinberg [16], Taylor window functions are from Collier 
{6], Balanis [18] and Pozar [19], and the Chebyshev Array is from Balanis [18] and Pozar 
[19]. 

The rectangular aperture usually abbreviated RECT has a uniform excitation over 
the array length L. This excitation is typically assumed to be unity. The magnitude of the 
directional gain pattern has a value of -13.5 dB at the first side lobe, the envelope halves 
with every doubling of spatial angle. 

The raised cosine and raised cosine squared taper functions have improved 
characteristics in that their first sidelobe levels are -23.5 dB for pedestal height of .5 and 


-32 dB for -edestal height .1 respectively. For a normalized array length p from -1.0 to 





1.0 and a pedestal height of a the element excitations are given by equations 2.18 and 
2.19. Note that for a selected pedestal height of 1.0 both these taper functions will reduce 


to a RECT window. 
1, = a+ (0.5 — 0/2) cos(pT]) (2.18) 
1,=a+(1—a)cos(pl) © (2.19) 
The Dolph-Chebyshev (Tschebysheff) array will create an array which has the 
narrowest main lobe in the look direction for a specified side lobe level. A unique 
characteristic of this array is that all the sidelobes are at the same level. The element 
excitations are related to Chebyshev polynomials. Several methods of efficiently 
computing the element weights of the Chebyshev array have been proposed. The method 
used in ARRAY, a nested product algorithm, proposed by Bresler [21] is recommended 
by Pozar [19] and Bums [20]. 


For a desired sidelobe level of Ry, in dB, R in volts is 


Rgy!20) 


R=10 


For an array of M elements the Chebyshev polynomial z, is given by equation 2.19. 


28 

















1 S 
Zy= cost sj cosh | (2.19) 


Defining the array element numbering scheme as depicted in Figure 2.7, N is the element 
number furthest from the array centerline. The array element weighting factors are 
calculated using equation 2.20. 


ee eae) for sie 2.20) 


1.0 for n=0 


p=[v-) for M 7 
~{N for M_ odd 


where: 


NP(n,a)= Lo "Tlf 
m=! 7 
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SN en ae 
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The Taylor weighted aperture, referred to as Taylor Line Source (one parameter), 
will produce a reduced sidelobe array which unlike a Chebyshev array does not have an 
optimum main lobe beamwidth to sidelobe level ratio. The array is very close to 
optimum though. It does require consideration as a design however, because unlike the 
Chebyshev array which has uniform sidelobe levels, the Taylor array has sidelobe levels 
which decrease monotonically as the angle increases from the broadside direction. This 
type of pattern might be advantageous when a known interfering source is located in a 
direction associated with the third or higher sidelobe. In this case the designer may be 
willing to trade some small amount of increased main lobe beamwidth for a lower 
sidelobe in the vicinity of the interfering source, a choice not available to a Chebyshev 


array. The Chebyshev array has sidelobes at a constant level. 
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For normalized array length p with values of -1 to 1, Taylor parameter n is a 
constant input by the designer which relates main lobe width to sidelobe level, and noting 
that is a modified Bessel function of order zero, the element excitations are calculated 


by equation 2.21. 


1,=1(nV1 -p’) | (2.21) 


Table 2.1, extracted from Balanis [18], documents the expected sidelobe level for the 
Taylor parameter n. 


Table 2.1 Taylor parameter (n) versus first sidelobe level. 


First Sidelobe Level (dB) 
4.0090 | 4.7551 | 5.4711 


The final window function incorporated into ARRAY is a user defined series of 






weights. This allows a user not familiar with beamforming theory to experiment with the 
weight values for each hydrophone position and graphically see the result. Additionally 
this selection allows weights from an externally calculated analysis to be incorporated 
into ARRAY where its graphics capabilities can be used to plot the resulting beam 
patterns. Possible candidates include Hann, Hamming, and Hansen-Woodyard End-Fire 
Array. 

Equations 2.14, 2.16, and 2.17 can now be used to analyze various array 
architectures including phase and amplitude tapering, by calculating and plotting the 
directional gain. These, along with the desired amplitude taper functions, calculated by 


equations 2.18 through 2.21 as applicable, form the main numerical analysis of the design 
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tool ARRAY discussed in detail in Chapter 4. ARRAY was used to conduct the initial 
study of array architecture within the constraints of the reverberant tank’s physical 
dimensions as well as the signal limitations imposed by the active control of scattering 


experiment. The output of the design study is included in Appendix B for review. 


2.3.4 Conventional Beamformer Output Power 
Often it is convenient to examine the average output power as a function of look 
direction @. For array output signal z(t), average power is defined as: 
P(®)=Elz(t)') 
The summation of the individual hydrophone signals, z(t) has been previously defined for 


i elements with weights w,(6), and received signal sp,(t) as: 
M 
z(t) = x w'(O)Sp,(t) 
Substituting and rearranging into a matr . formulation the average output power can be 
rewritten as: 
P(®) =wE[z(t)'z(t)]w 
For the conventional beamformer the weight applied to each element is the product of a 


phase taper to steer to a desired look direction, and an amplitude taper to control the 


sidelobe levels. 





2.4 Linear Array Synthesis 


Physical constraints imposed by the BBN acoustic tank as well as the array support 


structure limited the array aperture to approximately 10 feet. The active control of 
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scattering experiment had been designed to operate at a frequency of 2 kHz. These 
predetermined constraints allowed the array design to vary hydrophone separation, the 
positioning of the array , the amount of beam steering or phase taper, and finally the type 
of amplitude taper. In order to rapidly examine a large number of array architectures the 
design tool ARRAY was written and used extensively. Figures 2.9 through 2.12 are 


actual screen images captured while executing ARRAY during the initial design stage. 
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Figure 2.9 Eight element spacing analysis. 


Although the active pulse was predetermined to be a 20 msec sine wave at 2kHz it 
was decided to analyze the acoustic environment following scattering off the actively 
controlled plate up to 5 kHz. T.1is would allow the effect of the active control system to 
be better examined over a wider frequency range. At 2kHz, with sound speed of 4860 
feet per second (fps), wavelength (A) is 2.4 feet. Since spatial aliasing is not desired, an 
interelement distance of 1.2 feet is the maximum allowed. Keeping within the 10 foot 


aperture allowed results in an 8 element array. Figure 2.9 is the ARRAY screen output 
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showing the expected beam pattem for this first iteration at both 2 kHz and 5 kHz. For 
the initial iterations the beam was steered to broadside and a rectangular amplitude taper 
of unity applied. 

As is evident, interelement spacing of A/2 at 2 kHz results in dramatic spatial 
aliasing at 5 kHz and above. The solution requires an analog anti-aliasing low pass filter 
to eliminate any signal above 5 kHz prior to analog to digital conversion and an 
interelement spacing which precludes spatial aliasing at the highest expected frequency 
(5 kHz). Selecting 0.583 feet or 7 inches as an interelement distance, the 10 foot acoustic 
aperture will support a 15 element array. Figure 2.10 shows the expected beam pattem at 
the 2 kHz active pulse frequency as well as the highest frequency of interest, 5 kHz. 


These are considered optimum within the physical constraints imposed. 
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Figure 2.10 Fifteen element array analysis 


The design as proposed exhibits higher sidelobe levels than desired. Since the 


physical design limits preclude adding more elements, it is possible to reduce the sidelobe 
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levels by employing amplitude taper at the cost of widening the main beam. Figures 2.11 
and 2.12 demonstrate the effect of a Taylor window and a Chebyshev window on the 15 
element arrays of Figure 2.10. ARRAY was used to determine the element weighting 
factors for a Taylor parameter of 3 and a Chebyshev sidelobe level of 30 dB. Broadening 


of the main beam is clearly demonstrated. 


























ARRAY #1 ARRAY #1 BEAM PATTERN 
Er 20008 
beierity «rps: 4866 “ ® 
Sheeting Crit 8. es ¢  -1ef OS a A a 
iieve length cftr’ 2.43 v pep A ee Cees a 
Look Angie (deg): 98.6 . — 
Window’ CHEBYSHEY 
DirectivituycdB>: 10.848 - ~30 Sawa —— 
By beam wiaths a7. 8 + -48 1 1 TT fit 
6 pt iz 
é -sed If TT Ta Ti 
oe 
o i OR OR | OO A | | 
Freq vie ARRAY W2 2000 ARRAY #2 BEAM PATTERN 
elerity (fps>: 4860 











LF Se SS 6 cS SO 


Directivityt dB»: 9.4 
3 dB bean width: 1 








Gc Theta> «db» 


ITERATE 7? <y/n> 








16e 
Theta (degrees) 


Figure 2.11 Amplitude tapering at 2 kHz. 


The element weights calculated by ARRAY for creation of the beam pattems of 
Figure 2.11 at 2 kHz were manually entered into the 15 clement array at 5 KHZ to 
examine the array response at the higher frequency. This can be noted by the annotation 
of a USER defined taper window in the array table portion of tl e figure. 

Several other amplitude taper functions were analyzed and the summary of element 
weights is contained in Table 2.2 along with the value of the highest sidelobe in dB, the 
main lobe 3dB b:am width and the arrays calculated directivity index in decibels. 


Detailed output, including individual beampattern graphs, tabulated normalized 
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Figure 2.12 Amplitude taper at 5 kHz. 


directional gain values and amplitude taper weighting values for the arrays analyzed are 
included in Appendix B for review as desired. The data included in Appendix B was 


used to determine the effects on beam patterns due to changes in: 


Number of array elements. 
¢ Frequency response. 


* Steering the main response axis. 


Interelement distances. 


Amplitude taper functions. 

Appendix B contains example data from examining several iterations of array 
design. In actuality because of the ease with which a large number of beam pattems can 
be rapidly analyzed on screen many more types of array architectures were examined and 


discarded without producing hardcopy. 





Table 2.2 Amplitude window weights 


a 
a=.5 o=.1 n=3 -30 dB 
| 0.6302]. 


k : . 1 
: P ; 1 
1 
1 







fEenenss | _o6ss9 | 08306 | «ai00_| _oasar | 1.0000 | 
fEenew® [| 03750 _0s000 | 2aa04 | 05000 [1.0000 | 
fEiemen 9 | o7a52 | _o9ssa_| a7ooe | asvae | 10000 
feemen 10 | 06559 | oas06 | _4130_[ 08047 | 1.0000 | 


i 
siacive [2 | 23 | 26 | 0 | 





Based on 15 elements spaced 7 inches apart for 2 KHZ signal, c is 4860 ft/sec. 


Based on a review of this data the array employed consisted of 15 elements placed 7 
inches apart on an array support structure positioned to keep the broadside beam directed 
at the target plate, which eliminates the need to phase taper in order to steer the main 
beam to the desired target. The projector would be placed in the null between the first 
and second sidelobe, approximately 124 degrees off array axis. Post-processing the data 
would afford the capability of trying several amplitude weighting functions to determine 
which was best suited to the problem geometry, although a Chebyshev array appears 


optimum. Due to the fact that this array is specifically intended for spatial filtering of 
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multipath interfering signals and not extremely accurate direction finding, it is logical to 
use amplitude tapering and accept a wider main lobe in order to obtain the better 


reduction in sidelobe levels. 
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3 Experimental Setup 


3.1 BBN Reverberant Tank 


In water testing involved forming the synthetic aperture by successive hydrophone 
placement in the acoustic test facility at Bolt, Beranek, and Newman (BBN) located in 
Cambridge, Massachusetts. The tank facility is designed for measurement of acoustic 
power levels. The side walls and bottom of the tank were constructed so as to eliminate 
parallel opposing walls which could lead to standing wave formation. Figure 3.1 shows 


general physical dimensions of the facility. 


32 Fr 23 FT 
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Figure 3.1 BBN Reverberant Water Tank 


The reverberant tank holds approximately 11,000 ft’ of fresh water. Sound speed is 
estimated to be 4860 ft/sec based on density and temperature considerations. The tank 
bottom and four walls are lined with a highly reflective polyurethane foam which acts as 
a pressure release surface. The average absorption coefficient has been previously 
determined to be 0.05. The tank surfaces being highly reflective tend to create a large 
number of multipath signals, highly correlated and perhaps coherent, to the original 
projector pulse. These multipath signals are the predominate problem associated with 


beamforming and array processing in a reverberant tank. 
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3.2 Scattering Plate 


The target for the acoustic pulse is expected to be a steel plate with an attached 





control system which will cause the plate to respond as an extension of the free surface. 


Figure 3.2 shows the physical dimensions of the target plate. 
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Figure 3.2 Scattering Target Plate 


The plate is approximately 16 ft’ of 1/4 inch steel. The plate has applied damping 
treatments to establish a damping loss factor n of 1 to 2 percent. In order to allow the 
plate to be free floating, and to fluid load only one side, a cofferdam of 3/4 inch pine was 
built. The cofferdam was rigidly attached to the plate and sealed to prevent water entry. 


The plate, when waterborne, floated with a waterline of 3.5 inches as expected. 


3.3 Projector and Active Pulse 

The equipment used to generate the active pulse is depicted in Figure 3.3. The 
pulse consisted of a 40 cycles (20 msec) of 2 kHz sine wave with an amplitude of 4.5 
volts and an initial positive slope. A Hewelett Packard signal generator model HP 3314A 


was used to create the initial signal. This signal splits at the signal generator output with 
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one path being monitored on a dual trace oscilloscope and recorded as the trigger pulse 
on channei 3 of a TEAC XR-5000 14 channel tape deck. This signal will be the time 
synchronization for analog to digital conversion discussed under data collection and 


reduction. 
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Figure 3.3 Projector schematic 


The second signal from the signal generator passes through a single channel of a 
BGW 8000 power amplifier and a 19:1 transformer. The output of the transformer is 
monitored by a Beckman VOM to set the voltage to 150 volts using the power amplifier 
controls. The transformer output is sent to a Sumner & Mills model J-9 4 inch projector. 
The active pulse can be triggered manually or by time synch selected on the HP 3314. 
During the experiment the active pulse was automatically triggered at 10 second 


intervals. 
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3.4 Hydrophone Assembly 
The receiving section of the apparatus is shown in figure 3.4. One Celesco LC-10 
cylindrical hydrophone is positioned at 15 separate positions on the array support 


structure. 
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FROM PROJECTOR 
Figure 3.4 H-drophone schematic 


The hydrophone has an omnidirectional response orthogonal to its longitudinal axis 
and a cartoid response parallel to the longitudinal axis. The hydrophone has a response 
of -101 dB re lvolt/Pa. The signal is amplified by an ITHACO model 453 amplifier 
whose signal was monitored using headphones to detect any clipping. The amplifier was 
set to get an output amplitude close to + 5 volts as monitored on the dual trace 
oscilloscope. The signal was recorded on channel one of the TEAC XR-5000 
multi-channel tzpe recorder. Forty pulses were recorded for each cf fifteen hydrophone 
positions. Twenty were taken with tue iurget plate in the water and twenty with the plate 
removed from the tank facility. A continuous voice track was recorded to facilitate 
digitization and data transfer at MIT from the TEAC recorder to the MASSCOMP 


computer. 
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3.5 Forming a Synthetic Aperture by Hydrophone Position 

The acoustic aperture of the array was formed by sequentially positioning the same 
hydrophone at fifteen separate positions in the reverberant tank. In order to post process 
the data as if one 15 element array was used required reproduction of the active 
interrogation pulse and the capability to position the hydrophone accurately. The HP 
3314 signal generator was capable of creating repetitive pulses. In order to ensure 
accurate hydrophone positioning an array support structure was built. Figure 3.4 depicts 
the experiment as it was placed in the BBN tank. The projector, plate, hydrophone, and 


array support structure are shown. 
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Figure 3.5 Synthetic aperture array structure 


The array support structure is 20 feet of 2 inch radius PVC pipe free flooding with 
water. The pipe is divided into three sections, a 3 foot stub tube at the bottom, 10 feet of 
acoustic aperture, and 7 feet at the top for handling and support. On the 10 foot acoustic 


window portion, a 12 inch long, 2.5 inch radius movable collar is mounted. This collar 





can be positioned by downhaul and uphaul cables to position it along the 2 inch PVC 
pipe. To prevent azimuthal rotation of the movable collar it is slotted and moves along a 


track attached to the back of the support structure. A mechanical scale allows positioning 
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the collar to within 1/16 inch along the array support structure axis. One foot of 3/4 inch 
diameter freeflooding PVC pipe is rigidly mounted to the movable collar orthogonal to 
the array axis. The LC-10 hydrophone is attached to the end of this pipe to keep it away 
from the support structure. 

The entire array support structure is positioned in the tank so as to make a 60 degree 
angle between array axis and tank bottom. This "array steering" allows the broadside 
beam to be directed toward the target plate and eliminates the need to phase taper for 
beam steering. 

By simple geometry the plate and projector are positioned such that the signal 
arrives from 30 degrees below the horizontal axis and is reflected at the same angle in 
accordance with Snell’s law toward array window centerline. The plate center is 14 feet 
4 inches from the wall near the array. The projector is 10 feet 6 inches from plate center 
and 7 feet from the tank wall. The projector is placed at a depth of 6 feet 2 inches to 
establish an angle of 30 degrees with the target plate and to place the pre 2ctor at 124 
degrees referenced to the array longitudinal axis. Referring to Figure 2.9 and the 
tabulated values in Appendix B shows this should place the direct path signal from the 


projector in a null. 


3.6 Data Collection Reduction and Processing 

The hydrophone output data is collected at tre BBN acoustic tank and recorded on 
the TEAC XR-5000 tape deck. Figure 3.5 depi :ts the digitization of the data at MIT 
using the MASSCOMP computer. The taped hydrophone data is replayed through a 5 
kHz lowpass anti-aliasing filter. The AD12FA analog/digital (A/D) converter along with 
a SH16FA sample and hold module are used to digitize data. Two analog to digital 


channels are opened. Channel 0 contains the trigger signal originally recorded as output 
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from the signal generator. Channel one contains the hydrophone output. 
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Figure 3.6 Processing of taped data 

The data acquisition on the MASSCOMP computer takes place inside the LAB 
Workbench (LWB) environment. The analog input is demodulated to separate channel 0 
from channel 1. The channel 0 signal provides the input for a trigger module which will 
pass channel | (hydrophone) data only after the 1 volt rising trigger is detected for greater 
than | jtsec. The output of the trigger module is a time clipped section of digitized 
hydrophone data which is recorded in a user defined file. The digitized, clipped data is 
viewed on a digital scope display on the LWB screen. The time window selected runs 
from trigger detect to 60 msec. This is sufficient time to see the direct path signal as well 
as plate scattered signal and some interfering multipath signals. 

Each individual pulse is digitized, timeclipped and recorded in a file on the 
MASSCOMP. Figure 3.7 shows the virtual instrument created in the Lab Workbench 


environment used to digitize, demodulate, timeclip and record a single hydrophone pulse. 
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Figure 3.7 Lab Workbench A/D conversion virtual instrument. 


Appendix C1 contains a series of plots forming a library of the digitized 


hydrophone output signals recorded by this method. 
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4 Programs to Assist in Array Design 


4.1 Beamforming Analysis (ARRAY) 


ARRAY was written as a design tool to allow rapid iteration of array parameters 
and to present a graphical view of the resulting beam patterns. Two separate array 
architectures can be shown on the screen simultaneously to aid in this analysis. The 
program is intended to be used in the initial phase of array design where it is beneficial to 


examine many possibilities quickly so as to eliminate unsatisfactory designs before 





significant effort is expended in numerical analysis or array construction. ARRAY may 
also prove valuable for those cases where significant postprocessing of prerecorded data 
is done. ARRAY as currently written will determine the beampattern over 360 degrees in 
one plane for a linear array of equally spaced elements. The program incorporates phase 
taper to support beam steering as well as amplitude taper to control sidelobe levels. The 
code has the capability to be upgraded to plane arrays or non-uniform element spacing. 
The graphics routines used have the ability to plot 3-D beam patterns should this 
enhancement prove desirable in the future. 

The user has the option to enter the number of elements, signal frequency, 
interelement distance, sound speed, steer angle and amplitude taper window. The 
program has been compiled and run on a Video Graphics Adapter (VGA), an Enhanced 


Graphics Adapter (EGA) as well as Hercules Monochrome Graphics Adapter. 


4.1.1 Program Compilation and Use 


The source code for ARRAY is included in Appendix A. The program was written 
in Microsoft Quick C and compiled on an IBM PS/2 Model 70 386 personal computer. 


The code makes use of the Microsoft extensions to K&R C specifically with respect to 
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string handling and graphics, and therefore may require modification prior to being 
ported to non-DOS systems. The program was compiled using the medium memory 
model with a stack size of 16K bytes in release mode. 

The graphics routines used for plotting are modified Microsoft graphics drivers 
provided by Quinn-Curtis Software for Industry, Science and Engineering as source code. 
As such the source code is proprietary and not included in Appendix A. The final 
program ARRAY is based on linking five separate modules together, array.c, menu.c, 
beam.c, plot.c, and output.c. Array.c is the main calling program from which execution 
Starts. Menu.c is the prime executive which controls initialization as well as program 
execution. Calls to beam.c, plot.c and output.c are made from menu.c. Beam.c is the 
primary calculation module where the normalized array gains are calculated as well as 
the individual element weights based on a selected amplitude taper function. Plot.c 
controls the graphical user interface, establishing windows, plotting the array gains and 
updating the array data tables. Output.c opens output files for data recording and ports 
graphics to the printer as desired by the user. 

The program is executed by entering ARRAY at the DOS prompt. Internal to the 
main calling program ARRAY looks for the file array.con on the default drive. This file 
contains information regarding the type of graphics driver and monitor installed as well 
as the printer used. If array.con is found the program loads the hardware configuration 
data and shifts execution to the program executive, menu.c. Should array.con not be 
found as would be the case on the first use of the program, a warning is issued and the 
required configuration data is prompted via menu. The configuration file array.con is 
opened on the default drive, the data written and the file closed. Program execution 
continues as before. The file array.con is now available the next time ARRAY is used. 
To change the hardware configuration delete array.con prior to executing ARRAY. The 


Hercules version will first load the graphics driver MSHERC and then execute. The 
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initialization routine in Menu.c assigns values to both array #1 and array #2, calculates 
array gain and plots the graphs. Figure 4.1 is the first screen displayed and shows the 
initialization array values. Refer to Figures 2.8 through 2.11 for other examples of the 


actual screen displays while running ARRAY. 
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Figure 4.1 ARRAY initialization screen 


The command window in the lower left of the screen prompts as to the desirability 
of iteration. Entering a character (N) will exit the program. Entering a (¥) will cause the 
command window to prompt for which array to iterate. The user should enter the desired 
array number({1) or 2)) to be prompted for array parameters. The program terminates 
the graphics window format for rapid entry of array characteristics. The input screen 
shows the current default value of each parameter in brackets < > which can be selected 
by simply using (<-ENTER ) vice typing a value. Figure 4.2 shows the parameter entry 
screen following entry of the sound speed. 

During entry a menu of window choices is given. The user should select the table 


value (CL through (6)) of the desired choice. After entering sound speed, or accepting the 











default, if the window function is any but a rectangular window the user will be prompted 
for more information. Once the element amplitude weights have been calculated they are 
presented on screen for verification. The user has a choice of accepting these values or 


retuming to a prior level of data entry depending on the window. 
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Figure 4.2 Array data entry window. 

To accept the element weights as shown enter (Y) and the program returns to a 
graphics window environment showing the updated array normalized directional gain 
pattern as well as the new array data table now updated to show the revised array 
characteristics, amplitude taper selected and the computed values of mainlobe beamwidth 
and directivity index. Since the program displays two completely independent arrays the 
user can shift between each at random making iterative changes and examining the 
results. 

The command window again prompts for iteration and the program continues in 
this manner until no further iterations are desired. When {N) is entered to stop iteration of 
array designs the command window requests the users desire to produce output from the 


array synthesis. The user has the option to respond yes, (¥) or no, (N). Yes will cause the 


command window tc prompt for which array to use 1,2, or both. 











Output from ARRAY is directed to two separate devices. The printer declared in 
the configuration file array.con is used to produce a full page landscape orientation of a 


single array characteristic table and the associated graph of normalized array gain. 


Simultaneously, the program opens a file array.dat on the default drive and saves the 


array characteristics (ie frequency, number of elements etc...), the calculated element 
amplitude weighting factors, and the tabulated array gain G[6] at whole degrees between 
0° and 180°. This file is then available for future reference and/or printing as the user 
desires. 

After producing the desired output ARRAY returns to the windows environment. 
The command window questions if the user desires to quit. A no will return the program 
to array iteration. A yes will cause program execution to terminate and a retum to the 


operating system. 


4.1.2 Numerical Analysis 


Three numerical algorithms are employed in calculating the normalized array gain 
in ARRAY. If the window function chosen is the rectangular function it is numerically 
much more efficient to use equation 2.14 vice the summations in equations 2.16 or 2.17. 
Additionally, while not requiring recursive summation loops, equation 2.14 is applicable 
to either even or odd numbered arrays while equations 2.16 and 2.17 require determining 
if the array is even or odd numbered. However, if any window function other than a 
rectangular one is selected only equations 2.15 or 2.16 are applicable. ARRAY will 
chose equation 2.14 to calculate array gain if the RECT window is selected and either 


equation 2.16 or 2.17 as required for all other amplitude taper choices. 





Equations 2.18 through 2.21 are coded in beam.c to calculate the element window 
weights. For continuous functions the function value is sampled at the normalized 
aperture position, from -1.0 to 1.0. This value is assigned to the element weight. 

The Taylor window function discussed in Chapter 2 requires the value of a 
modified Bessel function of order zero. ARRAY incorporates the modified Bessel 
function algorithm discussed by Press [7]. The function is based on the coefficients of 
Abramowitz [22] and is included in beam.c as function bessi0( ) for review. 

The Chebyshev array element weighting factors are given by Bresler’s [21] nested 
product algorithm in equation 2.20. In order to obtain the correct weighting factors 
requires solving equation 2.19 for the Chebyshev polynomial z). The compiler used to 
create the executable version of ARRAY does not support hyperbolic and inverse 
hyperbolic functions. These functions were approximated by their exponential forms as 


found in Beyer [8]. Exponential and log functions are supported by Microsoft Quick C. 


cosh(x) =5(e" +e") : (4.') 
cosh ‘x = log (x +Vx?- 1) (4.2) 


Directivity and directivity index, discussed in Section 2.3.3, are one quantitative 
assessment of the effectiveness or an array. To calculate directivity requires integration 
of the normalized directional beam pattern over 360 degrees. ARRAY calculates array 
response in a plane vice over all physical space. Numerical integration takes place in 
beam.c during the call to function dir( ). The extended Simpson’s rule discussed by Press 
[7] was used to numerically integrate G[8] using the 720 calculated values over 360 


degrees. 


*N 1 4 2 4 2 4 1 1 
i Sf(x)dx nol fasgirdnesne. Gh ot$n.th}+o[ | (4.2) 
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N is the number of points used in the numerical technique (720), h is the spacing between 
these points (.5 degrees), f, the value of normalized directional gain pattern G[9] at point 
n, the limits of integration are from 0 to 360 degrees. This technique is a closed formula 
as the boundary values are evaluated. The formulation of Equation 4.3 in code results in 


summing overlapping groups of three values. Equation 4.4 demonstrates the technique. 
an Lec AA 
Hy f(x)dx = Aish tshtah 
1. 4, 1 
eat ae 
1, 4, 1 


This gives the required 2/3, 4/3 alternation throughout the interior of the summation 
equation while preserving the 1/3 factor on the boundary elements. The error expected in 
Equation 4.3 is on the order of 1/N’. 

Mainlobe beamwidth is quantified by the -3 dB beamwidth. The values of 
normalized directional gain are examined at 1/2 degree increments from the desired look 
direction (where G[8] = 1.0). At the point where G[6] is greater than -3dB the iteration 
ceases. A linear interpolation between the two points bounding the -3dB level is made. 
The actual beamwidth is twice this value as the beam extends in both directions away 
from the main response axis (MRA). While not exact, this method is accurate to within 
fractions of a degree and is certainly adequate for the type of initial design review for 
which the model was created. 

The normalized directional gain pattern is calculated from 0 to 180 degrees and 
recognizing the symmetry of a line array beampattern these values are simply reproduced 
at angles 181 to 359. The gain is calculated at half degree values . This appeared to give 


sufficient resolution to the output plots and improved the accuracy of the numerical 
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integration required for directivity calculations. Should a higher resolution be desired at 


a later time the angular step size could be decreased by a simple change to the source 
code: 

ARRAY uses dynamic memory allocation to create arrays and vectors of user 
defined length at the time of program execution. The normal method of using arrays and 
vectors in C is avoided as it would place unnecessary limitations on the program or use 
significantly greater memory than required for the majority of arrays analyzed. An 


excellent discussion of the rational for this choice is contained in Press [7]. 


4.1.3 Software verification 
Any computer model no matter how sophisticated looking is only as good as the 
results it calculates. As a check of the coding of the program the same problem could be 


executed using three separate paths through the code. A RECT window function which 





uses equation 2.14 should give the exact same results as a USER window function with 
element weights 1.0, which uses equation 2.16 or 2.17. Similarly, a COSINE on a 
pedestal window function with the fractional pedestal height & set to 1.0 returns a set of 
weights equal to the RECT function. Figure 4.2 is the ARRAY output which compares 
the RECT window with a USER defined window of magnitude 1.0. This exercise 
demonstrates that ARRAY using equation 2.14 or equation 2.16 will arrive at the same 
solution. The same comparison was made using COSINE window function with a 
pedestal height of 1.0. 

Function calls to bessi0() and dir() in beam.c calculate the modified Bessel function 
of order 0 and the integration of the curve of array factor by Simpson’s 3/8 method. Both 
these functions were verified by comparison to known solutions to simplified problems 


prior to incorporation in ARRAY. 
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Figure 4.3 Comparison of array gain formulation. 

In order to verify the performance of ARRAY, the program was used to verify 
solutions to array problems found in Collier [6], Pillai [13], Balanis [18], and Pozar [19]. 
The problems were chosen so as to check as many of the functions of ARRAY as 
possible. Table 4.2 is a list of some problems run as well as the source of the results used 


to validate ARRAY as a design tool. 


ARRAY was used to compute the amplitude taper element weights, the sidelobe 
level in dB, the directivity in dB, and the beamwidth of the main lobe. These values were 
compared to the source document analytical solutions where possible and the error noted. 
Sidelobe, beamwidth and directivity error are presented as absolute values. Because the 


element weights varied considerably the largest absolute error for any singie element was 


converted to a percent error. Table 4.3 presents this information. 

















Table 4.1 Array validation exercises 


PROBLEM STATEMENT SOURCE FUNCTION 
DOCUMENT{ __ VALIDATED 
PT] cere anekechaers | a) [Rec 
10 element, d=4/4, RECT, steer 60 ° | 18] |= RECT, Steer 


17 element, 


10 element, d=A/2, USER, steer 90 °, (Binomial element [18] USER window 
weights) 

10 element, d=/2, sidelobe level -26 dB, CHEBYSHEV, be CHEB YSHEV 
steer 90 ° 












19 element d=.7A, n=6, TAYLOR, steer 90 ° 


17 element, d=A/4, sidelobe level -30 dB, CHEB YSHEV, [18] CHEB YSHEV 
steer 90 ° 





~ 


Table 4.2 ARRAY output error 


Problem Weight Sidelobe Main lobe 
Error (%) Level error Beamwidth 
(dB) Error(*) 





*. 1.0 across window; NC- not calculated in source; NA-not applicable 


The numerical techniques used in ARRAY appear to be sufficiently accurate. The 
errors associated with amplitude taper weights, sidelobe level, and mainlobe beamwidth 


are extremely small. The errors noted in directivity are not unexpected. ARRAY 























numerically integrates the area under the normalized directional gain graph and ratios the 
value obtained to a unity gain level over 360°. The source documents rely on algorithms 
which calculate directivity for a broadside array only. ARRAY is more robust in that it 
will calculate directivity irrespective of the look direction. The difference in directivity is 
not large and the program demonstrates the correct trends in directivity (ie amplitude 


taper decreases the directivity index a smal] amount). 


5 Signal Processing 


5.1 Data Collation and Reduction 


Section 3.6 discussed the method of digitizing each pulse by using Lab Workbench 
on the MassComp computer. Figure 3.7 schematically showed the virtual instrument 
used. Each pulse of 60 msec duration was stored in an individual file. Each file 
contained 600 sample points which defined the pulse sufficiently to prevent aliasing up to 
5 kHz. For each of the 15 hydrophone positions which defined the synthetic aperture, 40 
pulses were recorded, 20 with the scattering plate in position and 20 with the plate 


removed from the tank. 
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Figure 5.1 Data reformat virtual instrument 


Each file was stored as a binary file on the MassComp. A second set of files stored 
in ASCII format were created to allow transfer from the UNIX based MassComp to a 
DOS based system. This allowed the raw data to be manipulated on a PC system and 
eased incorporation into documents. Figure 5.1 shows the virtual instrument used in Lab 
Workbench to reformat a data file from binary to ASCII prior to data transfer. 

The ASCII files created were ported to a NEC Powermate 386 PC based system via 


a serial port null modem connector using BITCOM communication software. The DOS 
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system files were then available for data manipulation and tabulation. Each set of 20 
pulses specific to a hydrophone position were averaged to obtain one representative 60 
msec pulse. This reduced 600 files, each with 600 values, to 30 files. 

Prior to actually carrying out the averaging manipulation the individual pulses were 
plotted on the same set of axis to verify each has a similar waveshape. Each file was 
available for amplitude weighting and summation as discussed in Section 2. The array 
was formed by phase and/or amplitude tapering and summation of these averaged files. 
The quantitative assessment of the array would be determined by a measurable difference 


between the array signal with the plate waterbome and that with the plate removed. 


5.2 Digitized Pulse Analysis 


Since it is easier to visually integrate plots rather than tables of data Appendix C 
contains the library of digitized pulse shapes examined throughout the experiment. Each 
graph plots the 600 points for each pulse as a function of time. The digitized output from 
the analog/digital converter uses a 12 bit word giving 4096 possible values (2'”). The 
analog input varies from -5 volts to +5 volts resulting in an A/D conversion of 0.0024 
volts per digitized unit value and values of digitized data from -2048 to +2048. The time 
scale is from 0 to .06 seconds or 60 msec. Figure 5.2 is a plot of the acoustic pulse at 
hydrophone 8 with the scattering plate in the water. As graphed the pulse is unaveraged 
and untapered. 

The 0 time point corresponds to the pulse trigger and the actual acoustic signal will 
arrive at the array some time later based on the direct path distance from projector to 
hydrophone. A review of Section 3 will verify that this distance is approximately 23 feet 
to the hydrophone at array centerline, hydrophone #8. Using 4860 feet/second as the 


sound speed in water, the acoustic pulse along the direct path should arrive .0047 seconds 
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Figure 5.2 Hydrophone 8 digitized pulse 


after trigger. Figure 5.2 supports this analysis. 

Appendix C1 contains a series of pulse graph. corresponding to unaveraged, 
untapered hydrophone data. For each hydrophone position the two cases corresponding 
to the scattering plate are presented. A cursory review shows that extraction of one 
specific signal path from the pulse is very difficult. Any difference due to the plate is 
marginal at best. Additionally, this series of waveshapes vividly demonstrates the impact 
of the reverberant tank. The projected signal is 20 msec of 2 kHz sine waves. The 
received signals at each hydrophone are all significantly corrupted by the phase addition 


and subtraction resulting from multipath signals. The number of, and effects due to these 


multipath signals are the most significant problem associated with conducting array 


processing in the BBN reverberant tank. 





A review of the 30 plots contained in Appendix C1 shows one potential problem. 
Hydrophone 11 is unique in that there is a significant difference between the case of 
scattering plate present and then removed. An initial evaluation might conclude that a 
multipath signal was responsible for almost perfect phase cancelation. It would be 
unlikely that the initial section of the direct path signal (4.7 - 5.7 msec) would also be 
canceled as it arrives at the hydrophone prior to any reflected path signal, an analysis 
tending to negate the proposed explanation. Additionally, it is expected that the free 
surface scattered signal should have a higher amplitude that the plate scattered signal. In 
order to examine the pulse, the scale of the axis were changed and the data replotted. 
Figure 5.3 shows hydrophone 11 with the scattering plate removed and the range of the 


digitized hydrophone output rescaled. 
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Figure 5.3 Hydrophone 11 digitized output 
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Comparing the pulse shape of Figure 5.3 to that of the same hydrophone with the 
scattering plate in the water contained in Appendix C1 shows that both have identical 
wave shapes. All twenty pulses used to average hydrophone 11 output with the scattering 
plate removed exhibit this phenomena. It is felt that during the data recording or 
translation phase of the experiment this signal was not amplified correctly. In order to 
minimize effect of the potentially corrupted data hydrophone 11 was averaged, tapered 
but not incorporated into the array beamforming. Since the major quantifiable 
assessment of the array is to examine the difference between array output with the 
scattering plate present and that with the plate removed, hydrophone 11 was not included 
in afray processing. 

Appendix C2 presents the same information in the same format but graphs the 
average of 20 pulses for each hydrophone. This series of plots when compared to 
Appendix C1 demonstrates the repeatability of each pulse. The pulse shape of the 
averaged signal is almost identical to that of each of the original pulses. As discussed in 
Section 3 the accuracy of a synthetic aperture is dependent on the repeatability of the 
active pulse. The averaging process is expected to minimize slight variations in the 
acoustic environment during the period of active interrogation and data recording. 


Appendix C3 contains the averaged hydrophone data from Appendix C2 tapered 


using the Chebyshev element weights calculated by ARRAY and presented in Appendix 


B as well as Table 2.2. This set of amplitude taper weights should result in -30 dB side 
lobe levels with a 17.9° main lobe 3 dB beamwidth for a signal at 2 kHz. The associated 


beampattern for any signal at 5 Khz is also included in Appendix B for review. 











5.3 Single Omnidirectional Hydrophone Analysis 


Section i discussed the possibility of using a single hydrophone appropriately 
positioned to get a quantitative assessment of the effectiveness of an active control | 
system. Because of the effects of the multipath signals the pulse length must be too short 
to give an adaptive control system sufficient time to respond. A 20 msec pulse while 
good for the control system results in a received signal like that in Figure 5.2. 

Simple geometric calculation will show that the reflected signal from the plate 
afrives at approximately .0050 seconds (varying slightly depending on hydrophone 
position) after trigger. Figure 5.2 shows that it is virtually impossible to separate the 
desired signal from the received pulse even with apriori knowledge of the arrival time. 
The scattered signal arrives at the hydrophone while the direct path pulse is still passing 
through. Compounding the problem, the rear wall reflection also overlaps the time 
window of the desired signal. Time clipping simply will not work when long pulses are 


required and signals which are reflections of the active pulse can interfere. 


5.4 Subtraction of the Reverberant Field 


A second method to quantify the scattered signal involves subtracting the 
reverberant field. To accomplish this the received signal of an omnidirectional 
hydrophone taken without the plate in the water is subtracted from the signal at the same 
hydrophone position when the plate was present. The assumption is that the difference 


between the received pulses is due to the effect of the scattering plate. 
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Appendix D1 contains a series of pulse graphs based on this method. Each plot is 
the difference between the averaged received signal for a specific hydrophone position 
with the plate present and with the plate removed. As discussed above the scattered 
signal of interest is expected to arrive at the hydrophones in the vicinity of .0047 to .0053 
seconds depending on hydrophone position. 

A review of the plots contained in Appendix D1 indicates that this method does not 
result in an easily discemable pulse attributable to the scattered signal (5 to 25 msec 
along the time axis). The significant reduction in the value of the digitized hydrophone 
output indicates that the signal difference between the two cases is small. For long pulses 
in a reverberant facility the method of subtracting the reverberant field does not seem to 
provide a difference in received signals usable to identify the effects of an active control 


system. 


5.5 Beamforming Witi Untupered, Averaged Signal 


The beamforming process uses equations 2.9 and 2.15 to phase and amplitude taper 
each hydrophone signal and then sum the results. Since the array was physically 
positioned to direct the broadside beam to the target plate phase taper is not required. 
Using a rectangular amplitude weighting function of unity should result in the normalized 
array gain pattern of Figure 2.9 with the array characteristics in Table 2.2. This 
information is also contained in Appendix B. The array processing should result in 
coherent addition of the scattered signal and incoherent addition of noise and signals 
from the non-look directions. The expected first sidelobe should be reduced by -13 dB 
and the -3 dB beamwidth along the main response axis is 14". 

Figure 5.4 shows the resultant array output signal after processing with the 


scattering target plate in the water. There is reinforcement of the signal in the 5 msec to 
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Figure 5.4 Untapered array output 


30 msec range as was desired. Additionally as array theory indicates, there is an 
amplitude reduction relative to the desired signal outside of the time window established 
for the arrival of the scattered pulse. The initiation of the pulse on the graph occurs at 5.5 
msec which is consistent with the geometry of the situation for a reflected path and later 
than the 4.7 msec pulse initiation seen in Figure 5.2 which is consistent with direct path 
arrival. There is some noticeable modulation of the pulse in the 5-25 msec window 
which is likely due to interference from alternate path signals. Considering the cone 
shape of the passive beam formed by a linear array this effect is expected. 

Appendix D2 contains pulse plots of the processed array output wi. tne plate both 
in and out of the tank. The pulse shapes are similar and show a slightly stronger signal 
when the scattering plate is absent. The results indicate that even simple non-amplitude 


tapered beamforming of a judiciously timeclipped signal will produce an output usable to 




















quantify the effectiveness of an active scattering control system on the target plate. The 
difference in magnitude of the pulse amplitude is difficult to see unless the pulses are 


plotted simultaneously. 


5.6 Chebyshev Array 


7500 


DIGITIZED HYDROPHONE OUTPUT 








TIME (SECONDS) 


Figure 5.5 Chebyshev tapered array output 


Array processing was repeated using the Chebyshev array amplitude weights of 
Table 2.2. These weighting factors should produce the beampattern shown in Appendix 
B. Figure 5.5 reproduces the plot of the tapered array output. The sidelobe levels all will 
be -30 dB with a -3 dB beamwidth of 17.9° on the main response axis. The results of this 
analysis are comparable to that of the untapered array and support the ability of array 


processing to allow the use of longer pulses in the active control of scattering 
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experiments. 


Appendix D3 contains the plots of the processed array signal with the plate in and 


out of the reverberant facility for comparison. The pulse shapes and indicate phase 





addition at the desired time (5 to 25 msec) and phase cancellation of signals not from the 
desired look direction. The pulse begizis to rise at 5.4 msec as expected for the reflected 
signal. The same modulation of the signal in the 5 to 25 msec window as was noted in 
the previous section is repeated. The amplitude of the signai when the plate is absent is 
noticeably larger than in the previous case. This is most likely due to the better sidelobe 
suppression expected in a Chebyshev array. 


5.7 Power Spectral Density 


In order to better visualize the effectiveness of an active scattering controller 
attached to the target plate it may be advantageous to examine the differences in power 
spectral density of the array output signal for several cases of interest. Such an analysis 
may result in a better interpretation of the effect the active controller has in a frequency 
band. As a validation of the array processing the power spectral density with the plate as 
a target can be compared with the same calculation with the plate absent. This would 
give an indication of the maximum change in the scattered signal that could occur should 
the active control system make the plate appear as a continuation of the free surface. The 
Lab Workbench environment can be used to calculate the power spectral density of the 
array output and to display the result graphically. 

A singie omnidirectional hydrophone, the untapered averaged array output and the 
Chebyshev array output were ported back to the MassComp computer by reversing the 
process discussed in Section 5.1. In the Lab Workbench environment the power spectral 
density was calculated, displayed and recorded. The virtual instrument created to 


accomplish this is schematically shown in Figure 5.6. 
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Figure 5.6 Power Spectral Density virtual instrument 


The resulting plots of power spectral density are included in Appendix D4. The 
plots include both conditions of scattering plate position with a single hydrophone, the 
RECT window array and Chebyshev tapered array architectures. The difference in power 
spectral density for the two cases supports the contention that array beamforming will 
differentiate plate scattering from free surface scattering with pulse lengths as long as 20 
msec. As expected the Chebyshev array is slightly superior to the rect weighted window. 
Table 5.1 summarizes the values extracted from the power spectral density graphs in 


Appendix D4. 
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Table 5.1 Power spectral density summary 


Processing Method Maximum Amplitude Frequency at Maximum 
(Volts’/Hz) Amplitude Hz 


Omni hydrophone w/plate 0.16296 
Omni hydrophone w/o plate 0.17670 2011 


RECT window array w/ plate 2.25710 1982 
i 2.81620 


RECT window array w/o plate 1982 
Chebyshev array w/ plate 1.12230 1982 
Chebyshev array w/o plate 1.44990 


5.8 Summary 





The ability to determine the change in the acoustic signal scattered off a plate 
located at the free surface is crucial to the study of the active control of scattering. 
Logistically and economically it is advantageous to conduct these experiments in a 
reverberant, acoustic tank. Due to the large number of multipath and reflected signals, 
the received pulse at one omnidirectional hydrophone can be significantly altered. A 
review of the received waveform at fifteen omnidirectional hydrophones documents the 
inability to accurately extract the desired signal scattered off the target plate. An attempt 
to simply subtract the reverberant field leaving only that part of the signal due to the plate 
was also ineffective. 

The formation of a fifteen element array, properly phase tapered to position the 
main response axis on the target plate was successful in differentiating between the plate 
in the tank and the plate removed from the tank. The use of an interactive beamforming 


model, ARRAY, was integral in determining array characteristics and amplitude 
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weighting factors used in the post processing of the recorded data. The standard 
untapered array with -13 dB main sidelobe was successfully demonstrated as was a 
Chebyshev array with -30 dB sidelobes. 

Creating the array aperture synthetically, by sequentially positioning the same 
hydrophone in fifteen different but known positions, was successfully demonstrated as a 
viable alternative to the costly method of building an array using fifteen separate 
hydrophones. Building a multi-element array would require testing to determine accurate 


phase and amplitude matching of the hydrophones prior to use. 
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6 Future Designs for Increased Capability 


This report documents the successful use of a linear array to detect changes in 
scattered acoustic energy from a plate located at the free surface in a reverberant tank. 
As the desire to conduct acoustic testing in the reverberant tank increases so will the need 
to use more advanced means of array processing. If economically feasible, the next step 
should be to eliminate the synthetic aperture and build multi-element arrays which could 
be processed simultaneously. The MASSCOMP computer and TEAC data recorder have 
the capability to conduct this type of experiment once the required hydrophones and 
amplifiers are purchased. 

In terms of array types the linear array should be replaced by the two dimensional 
plane array. This can significantly improve spatial filtering. The design tool ARRAY has 
the capability to be modified to allow plane array beamforming even to the extent of 
using different amplitude weighting functions in each dimension. The plotting routines 
are capable of contour and three dimensional plotting. Plane arrays will improve the 
signal to noise ratio allowing smaller scattering target fixtures. 

The plate used as a scattering target in this study was a 4 ft X 4 ft steel plate. The 
plate anticipated for use in the active control of scattering experiment may be as small as 
10 inches X 10 inches. This presents a significantly smaller target. The ease with which 
this study traded main response axis beamwidth for reduced sidelobe levels may no 
longer be a prudent consideration. The small effective cross section will require efforts 
to minimize beamwidth. Plane arrays and linear arrays of more elements may be 
required to obtain results comparable to this study. However, the physical dimensions of 


the BBN tank place an upper limit on how large the acoustic aperture can be. 
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High resolution beamforming methods should be investigated. Techniques such as 


Capon’s (9], maximum likelihood method may provide improverrents over classical 
beamforming techniques which justify their more sophisticated numerical techniques. 
The associated but separate techniques of adaptive signal processing discussed by 
Widrow [14] may be useful in the environment of the experimental tank facility, where 


the interfering signals are repetitive and essentially unchanging. 
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Appendix A Interactive Beamform Program ARRAY 


74 


/* ARRAY.C,; MASTER CALLING PROGRAM TO CREATE AN ARRAY SYNTHESIS PROGRAM 
USED TO CALCULATE AND GRAPHICALLY PRESENT PLOTS AND TABLES OF NORMALIZED 
DIRECTIONAL GAIN PATTERNS FOR DISCRETE LINEAR ARRAYS OF OMNIDIRECTIONAL 
ELEMENTS WITH USER SELECTED WINDOWING FUNCTIONS. */ 

/* rev 1:2 APR 90 add TAYLOR WEIGHTING */ 

/* rev 2:5 APR 90 add CHEBYSHEV WEIGHTING */ 

/* rev 3: 12 APR 90 add directivity and beamwidth */ 

/* rev 4: 17 APR 90 include switch for various video adapters */ 

/* rev 5: 18 APR 90 add output to file/printer support */ 


/* PREPROCESSOR DIRECTIVES */ 
#include <stdio.h> 

#include <stdlib.h> 

/* GLOBAL VARIABLE DEFINITIONS */ 
/* EXTERNAL VARIABLES DEFINED */ 


int mi,m2: /* Number of elements in Array 1, 2 */ 
float f1,f2; /* Signal frequency in Hz areray 1,2 */ 
float d1,d2; /* Interelement spacing in feet array 1,2 */ 
float thetal ,theta2; /* steer “look” angle broadside = 90 */ 
float lambda I ,lambda2; /* signal wavelength */ 

float c; /* sound speed in fps */ 

int wint! ,wint?2; /* window type */ 

float *G1,*G2; /* pointer for unit offset vector gain */ 
float *winl,*win2; /* pointer for vector of f amplitude tapers */ 
int i, j, k: /* loop counters 

float dil ,di2; /* directivity for arrays 1, 2 */ 

float bw1, bw2: /* 3 dB main lobe beamwidth */ 

float vmode,pm; /* video mode and printer config flags */ 
char inline[80]: /* used for config data entry */ 

int nent; 


/* FUNCTION DEFINITIONS */ 

void menu(); 
/* MAIN PROGRAM LISTING */ 
void main() 

/* initiate main program */ 

OBEROI BOB IO ISI IO OIG IS III IE I II OISIOI SIO IOIDIDI SIDI OI IOI IOC Itt 
* OPEN FILE TO OBTAIN CONFIGURATION DATA ON VIDEO SYSTEM AND i 
* PRINTER HOOKUP 
nbn bensannsay end anuusnsssdoebesnecedenbebcansaceiwssshssessessnsennesee4) 
FILE * fp: 
FILE *out; /* ot data file pointer */ 
fp = fopen("array.con’,"r") 
if(fp==NULL) 

{ 


system("cls"); 
printf("*** ARRAY HAS NOT BEEN CONFIGURED FOR SYSTFM HARDWARE ***\n"); 
printf("*** FILE ARRAY.CON NOT FOUND IN DEFAULT DIRECTORY ***\n"); 
fp=fopen("array.con”,"a"); 
mode:printf(\nSELECT AND ENTER GRAPHICS ADAPTER FROM MENU BELOW:\\n"); 
printf(’(N1- CGA 4 Color\n\2- EGA\nN3- VGA\n\4- HERCULES\n"); 
printf('‘WNENTER GRAPHICS ADAPTER :"); 
gets(inline); 
Stofa(inline, & vmode,&ncnt, 1); 
switch((int) vmode) 
{ 
case |: 
vmode=1.0; 
break: 
case 2: 
vmode=6.0: 
break: 














case 3: 
vmode=8.0; 
break; 
case 4: 
vmode=12.0: 
break; 
default: 
printf(" INVALID CHOICE; PICK 1, 2, 3, OR 4Nn"); 
goto mode; 


system("cls"); 
printf(‘\nSELECT AND ENTER PRINTER FROM MENU BELOW\\n"); 
printf("\1- EPSON MX\nN2- EPSON LQ\nN3- TOSHIBA 24 PIN\o\4- HP LASER JET\o\5- HP THINK 
JET\N6- EPSON FX\n"); 
gets(inline); 
stofa(inline,&pm,&ncnt,1); 
m= P -1.0; /* shift indices to agree with InitSeGraphics() */ 
rintf(fp,"%f Jof",vmode,pm); 
fclose(fp); 


] 
if(fp!=NULL) 
{ 
fscanf(fp,”"%f %f",&vmode,&pm); 
} 


POI BAAIIIGI II IOI IOI IOI k 


* CONFIGURATION COMPLETE SHIFT TO PROGRAM EXECUTIVE MENU ( ) 4 


BEERS IO OCGA IO OI TO RA TTT 
menu(); 


/*end main program = */ 


/* MENU1.C LISTING TO CREATE A MASTER EXECUTIVE ENVIRONMENT FOR ARRAY.C */ 
/* PREPROCESSOR DIRECTIVES */ 

#Hinclude <stdio.h> 

#include <math.h> 

#include <stdlib.h> 

#include "worlddr.h" 

#include "segraph.h" 


#define PI 3.1415926 
/* EXTERNAL VARIABLES DEFINED */ 


int m1,m2; /* Number of elements in Array 1, 2 */ 
float f1,f2; /* Signal frequency in Hz areray 1,2 */ 
float d1,d2; /* Interelement spacing in feet array 1,2 */ 
float theta] theta2; /* steer "look" angle broadside = 90 */ 
float lambdal ,lambda2; /* signal wavelength */ 

float c; /* sound speed in fps */ 

int wintl ,wint2; /* window type */ 

float *G1, *G2; /* pointer for unit offset vector gain */ 
float *win!,*win2; /* pointer for vector of amplitude tapers */ 
float fm1, fm2, fwintl,fwint2; /* floats of integer variables */ 

int i, j, k; /* loop counters */ 

int init; /* flag for initialization of plot windows */ 

float dil,di2; 

float bw 1,bw2; 

float vmode, pm; /* video and printer configuration */ 
/* FUNCTION PROTOTYPES */ 

void plot(): 

void initial(); 


void menu(): 

float *vector(): 
void nrerror(): 
void free_vector(); 


void menu() /* MAIN EXECUTIVE FOR ARRAY SYNTHESIS */ 


{ 

char choice, chstr{2] ={0 ,0}, inline(81]; 

int nent, mode,num; 

/* ALLOCATE MEMORY FOR UNIT OFFSET VECTORS */ 
Gl=vector(0,720); /* DIRECTIONAL GAIN OVER 360 DEGREES */ 
G2=vector(0,720): /* DIRECTIONAL GAIN OVER 360 DEGREES */ 


/* Initialize arrays and compute directional gain factors 
and plot to screen */ 
initial(): 
/* OPEN COMMAND WINDOW, INQUIRE IF ITERATION DESIRED */ 
do /* master loop ends with quit program */ 
{ 
do 
{ /* open main array iteration loop */ 
SetCurrentWindow(1); 
ClearWindow(): 
if(vmode!=12.) 
SetPlotBackground(2); 
BorderCurrent Window(3); 
SelectColor( 15); 
settextstyleXX(0.0,1): 
settextjustify XX(0,0): 
Title Window("COMMAND WINDOW"), 
setlinestyleXX(0,0,3): 
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LabelGraphWindow(20,300,"ITERATE ARRAY <!/2>",0,0); 
choice = getch(); 

chstr[0]=choice; 

LabelGraph Window(850,300,chstr,0,0); 


CloseSEGraphics(): 
| haclihasiaaindiaeaiadiadiadbaabaediadiadlaaiatiadiaciediaaaiaeiiadladhadeaiediadladiacedtaclacladia diac dia ciate dladiatia ecient alecledladiaciadiedlecladladealediadlasiadiadisliclasiadiad 
* DATA ENTRY SECTION TO ITERATE ARRAY PARAMETERS * 


AAI CIOIOI ICICI OI IOI SIDI I III ISIS SII I IOI SI IRI SIDI IOI III IIB IIDIO III DIO I DITIO AT IA A IA | 
system("cls"); 
printf("\nENT ER THE REQUIRED INFORMATION HIT RETURN TO ACCEPT <DEFAULT>\n\n"); 
if(choice == ’1’) 


{ 
printf(/“ARRAY NUMBER ONE CHARACTERISTICS\n\n"); 
printf(‘(WN NUMBER OF ELEMENTS <%d>: ",m1); 
gets(inline); 
Stofa(inline,&fm1,&ncnt,1); 
m1=(int) ceil(fm 1); 
winl=vector(1,m1); 
~ printf(‘WN SIGNAL FREQUENCY (HZ) <%6.0f>: ",f1); 
gets(inline); 
Stofa(inline,&f1 ,&ncnt, 1); 
printf('W\t INTERELEMENT SPACING (FT) <%5.2f>: ",d1); 
gets(inline): 
Stofa(inline,&d1,&ncnt, 1); 
printf(‘WW STEER DIRECTION (DEG) <%4.1f>: ",thetal); 
gets(inline); 
Stofa(inline, &thetal .&ncnt.1): 
printf(’\N INDICATE DESIRED AMPLITUDE TAPER FUNCTION \n"); 
printf(’'\***CURRENT PROGRAM REQUIRES SYMMETRY ABOUT ARRAY CENTER***n"); 
printf(’/NNN I- RECTANGULAR‘); 
printf(‘\W\N 2- USER DEFINED\n"), 
printf("NNWN 3- RAISED COSINE\n"); 
printf(/NNN 4- RAISED COSINE SQUARED\n"); 
printf(/WNN 5- TAYLOR\n"); 
printf(/NNN 6- CHEB YSHEV\n"); 
printf(‘W““ ENTER MENU NUMBER FOR TAPER FUNCTION <%d> ",wint1); 
gets(inline); 
Stofa(inline, &fwint! ,&ncnt,1); 
wintl=(int) fwint!; 
printf('\\Wt ENTER SOUND SPEED (c) <%5.1f>: ",c); 
gets(inline); 
Stofa(inline,&c,&ncnt, 1); 
lambdal = c/f1; 


} 
if(choice == ’2’) 


{ 

printf(“ARRAY NUMBER TWO CHARACTERISTICS\n\n"); 
printf('\\t NUMBER OF ELEMENTS <%d>: ",m2); 
gets(inline); 

Stofa(inline,&fm2,&ncnt,1): 

m2=(int) ceil(fm2); 

win2=vector(1,m2); 

printf('(NN SIGNAL FREQUENCY (HZ) <%6.0f>: ",f2); 
gets(inline): 

Stofa(inline, & f2,&ncnt.1); 

printf(\WN INTERELEMENT SPACING (FT) <%5.2f>: ",d2); 
gets(inline); 

Stofa(inline,&d2,&ncnt,1); 

printf(‘\N STEER DIRECTION (DEG) <%4. 1f>: ",theta2); 
gets(inline); 

Stofa(inline, &theta2,&ncnt, 1); 

printf(NN INDICATE DESIRED AMPLITUDE TAPER FUNCTION \n"); 
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printf("(\***CURRENT PROGRAM REQUIRES SYMMETRY ABOUT ARRAY CENTER ***\n"): 
printf(/(NNN. 1- RECTANGULAR\o"); 

printf(‘NNNt 2- USER DEFINED\n"); 

printf('(NNN 3- RAISED COSINE\n"); 

printf(‘NNN 4- RAISED COSINE SQUARED\n"); 

printf(‘WNN 5- TAYLOR\no"); 

printf(‘NNN 6- CHEB YSHEV\n"); 

printf‘ ENTER MENU NUMBER FOR TAPER FUNCTION <%d> ",wint2); 
gets(inline); 

Stofa(inline,&fwint2,&ncnt, 1); 

wint2=(int) fwint2; 

printf(’NN. ENTER SOUND SPEED (c) <%5.1f>: ",c); 

gets(inline); 

Stofa(inline,&c,&ncnt, 1): 

lambda2 = c/f2; 


POR ORI IG IOI OI I I io IO i gO TOO IOI iio IOI Oi i foioi doioi ie io te 
* INITIATE GRAPHICS MODE AND PLOT REVISED ARRAY BEAM PATTERN * 
BEEBE RICE IO ISIE EIS CICISICE ICICI IIIS S IIE OI IID OR ii ttt i tet iio / 
if(choice == 1’) 
beam(m1,d1 lambda! theta! ,wintl,G1,winl &dil,&bw1); 
else 
beam(m2,d2.lambda2,theta2,wint2,G2,win2,&di2,&bw2); 
InitSEGraphics((int) vmode): 
plot(2,1,wint!,G1); 
plot(2,2,wint2,G2): 
free_vector(win1,0,m1): 
free_vector(win2,0,m2): 


SetCurrentWindow(1); 
ClearWindow(); 
if(vmode!=12.) 

SetPlotBackground(2); 
BorderCurrent Window(3); 
SelectColor( 15); 
settextstyleXX(0,0,1); 
settextjustify XX(0,0): 
TitleWindow("COMMAND WINDOW"); 
setlinestyleXX(0,0.3); 
LabelGraphWindow(20,600,"ITERATE ? <y/n>",0,0): 
choice = getch(); 
chstr[0]=choice: 
LabelGraphWindow(850,600,chstr,0,0); 
} while (choice != ’n’): /* close array iteration loop */ 
Clear Window(); 
LabelGraphWindow(20,600,"OUTPUT DATA ? <y/n>",0,0); 
choice = getch(); 
chstr[0]=choice; 
LabelGraphWindow(850,600,chstr,0,0); 
if(choice==’y’ll choice==’ Y’) 


{ 

ClearWindow(): 

LabelGraphWindow(20,600,"Array 1 2 both? <1/2/b>",0,0); 
choice = getch(); 

chstr[0]=choice: 

LabelGraphWindow(850,400,chstr,0,0); 

switch(choice) 


case ’1’: 


output(fl .c.m1,dI lambda theta! ,wint! dil ,bw1,winl,Gl); 
break: 
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case ’2’: 
output(f2,c,m2,d2,lambda2,theta2,wint2,di2,bw2,win2,G2); 
break: 

case ’b’: 
output(fl,c,m1dJ lambda! thetal wint! dil ,bwi,win!,G1); 
output(f2,c,m2,d2,lambda2,theta2,wint2,di2,bw2,win2,G2); 
break: 


} 
plot(2,1,wint!,G1); 
plot(2,2,wint2,G2): 


} 
SetCurrentWindow(1); 
ClearWindow(); 
if(vmode!=12.) 
SetPlotBack ground(2); 
BorderCurrent Window(3); 
SelectColor( 15); 
settextstyleXX(0,0, 1); 
settextjustify XX(0,0): 
TitleWindow(""COMMAND WINDOW"); 
setlinestyleXX(0,0,3); 
LabelGraph Window(20,600,"QUIT ? <y/n>",0,0); 
choice = getch(); 
chstr(0J=choice; 
LabelGraphWindow(850,600,chstr,0,0): 
} while (choice != ’y’); /* close array master loop */ 


CloseSEGraphics(); 
free_vector(G} ,0,720): 
free_vector(G2,0,720); 
/* end menu() the main executive program */ 


PAGO ICIOI BIC IOI RO ol i i UIOi iO IOI SOI ioioi ga IOS ICO IOI OI IOI aio i i i I IOI IOI ie te ieee i teak 


* CALCULATE ARRAY PARAMETERS AND PLOT DIRECTIONAL GAIN USING 
* INITIALIZATION VALUES FOR ARRAYS 1 AND 2, (see ARRAY.C) 


Ht eget NT SR RHHGRRAAAGSRLLAEHEISESERLAASRNAESTOREGASENEROESEHAAORANE«®/ 
void initial() 


im mode; /* set screen mode to single or dual plot */ 

int num; /* set array to be plotted */ 

/* Initialize both arrays: 2KHZ, 10/20 elements, steer 90/45 degrees, 1/2 
wavelength interelement distance, celerity 4860 fps, with rect window 
function of amplitude 1. */ 

ml=10; 

fm 1=10.; 

m2=20: 

fm2=20.; 

c=4860.: 

f1=2000.: 

f2=2000.: 

lambda! =c/f1; 

lambda2=c/f2; 

di=lambdal/2.: 

d2=lambda2/2.: 

thetal=90.; 

theta2=45.: 

wintl =]; 

wint2=]; 

fwintl=L.: 

fwint2=1.; 
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mode=2; /* set screen to dual plot mode */ 
if(vmode ==12.000) 

system("msherc’”’); 
printf(‘Initiating Default Graphics Window’): 
InitSEGraphics((int) vmode); 
beam(m1,d1 lambda] ,thetal ,1,G1,win!,&dil,&bw1); 
plot(mode,1,wint!,G1); 
beam(m2.d2,lambda2,theta2,1,G2,win2,&di2,&bw2); 
plot(mode,2,wint2,G2); 


) 
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| ha ddashdashital dadadtdadedhatsededeaededndadedidadedededesndedededndndnsnindnindadnaaaadeiadiianananindadn dla ddiabalahdadiebstahehilcesaaalel 


* BEAM.C CALCULATES THE NORMALIZED DIRECTIONAL GAIN PATTERN 
* G[THETA] FOR ANGLES 0<theta<179 CONSIDERING AMPLITUDE TAPER 
* THEN COMPLETES VALUES FOR ANGULAR SPACE BASED ON SYMMETRY. 


He ee Re ee ee ee Ae Ae Be A a A AR RR A RE EERE EE ER RAE EEE AERA EERE E EEE EE | 


/* PREPROCESSOR DIRECTIVES */ 
#include <stdio.h> 

#include <setjmp.h> 

#include <math.h> 

#include <stdlib.h> 

#define PI 3.1415926 

/* EXTERNAL VARIABLES DEFINED */ 


int m. /* Number of elements in Array 1, 2 */ 
float f: /* Signal frequency in Hz areray 1,2 */ 
float d: /* Interelement spacing in feet array 1,2 */ 
float theta: /* steer "look" angle broadside = 90 */ 
float lambda: /* signal wavelength */ 

float c: /* sound speed in fps */ 

int wint: /* window type */ 

float *G; /* pointer for unit offset vector gain */ 
float *win; /* pointer for vector of amplitude tapers */ 
int i, j. k; /* loop counters */ 

int wint; /* amplitude taper type marker */ 
char inline{8 1}; /* data ertry */ 

char ncnt: 

float *di: 

float *bw; 
/* FUNCTION PROTOTYPES */ 

void beam(); 

float bessi0(): 


float *vector(): 
double *dvector(): 
void nrerror(): 
void free_vector(): 


float dir(): 

float dbdown(): 

void beam(1n,d,lambda,theta, wint,G,win,di,bw) 
int m,wint: 

float d.lambda,theta,G[], win[],*di,*bw; 

{ /* open beam() function */ 
char ch: 


float num, den, ftheta{721], norm, fact, weight: 
float height. proj, x; /* window height and aperature projection */ 
double cheb.cheb!, side, angle: /* chebyshev polynomial, sidelobe */ 
int flag, thetai: 
double *a; /* chebyshev temporary weights */ 
thetai= (int) ceil(theta): 
theta = theta*PI/180.; 
if(wint==1) —/* loop for rect window uses [sin mx/(m sinx)] for speed */ 
{ 
norm=0: 
for (i=0;i<=360:i++) /* for spacial angles 0<theta<179 */ 


{ 
angle= i*PI/360.; 
den =PI*d/lambda*(cos(angle)-cos(theta)): 
ftheta[iJ=1.: 
if( fabs(den)>.00001) 
{ 
ftheta[i]J=sin(m *den)/(m*sin(den)); 


} 
ftheta[i]=fabs(ftheta[i]): 
if(ftheta{i]>norm) 
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normm=ftheta{i J: 
} 
} 


for(i=36 1.i<=720:1++4) 


{ 

ftheta[i]=ftheta[720-i]; 
for(i=0;i<=720:1++) 

{ 

ftheta[i] = ftheta[i]/norm; 

ftheta{ij=fthetali]* ftheta[i}: 

if( ftheta[i]<.0000001 ) 


ftheta[i] = 0000001. 





} 
G[iJ=10*log10(ftheta[i]): 


/* Determine 3 db main beamwidth & DI */ 
*di=dir( ftheta): 
*bw=dbdown(G.theta.’: 


} /* close conditional for rect window function */ 
if(wint == 2) /* user defined amplitude weighting window */ 
{ 
system("cls"): 
printf WARRAY AMPLITUDE WEIGHTS\n\n’): 
printf’ WWW NUMBER OF ELEMENTS <%d>\n "\m); 
repeat: for(j=1.j<=m:j++) 
{ 
printf('\nw ELEMENT %d WEIGHT : "j): 
gets(inline): 
Stofa(inline.& win[j].&ncnt, 1): 
} 


pontf( AWN VERIFY ENTERED WEIGHT VALUES\n\n"); 
forj=1:j<=m:j++) 


printf('\ ELEMENT %d WEIGHT : %8.4f\n ",j,win[j]); 


) 
pnotf("ARE THESE VALUES CORRECT <y/n>? "): 
ifich = getche() == 'n’) 


goto repeat: 


} /* close user defined amplitude window beam calculations */ 
if(wint == 3 Il wint == 4) /* RAISED COSINE AMPLITUDE WEIGHT CALCULATION */ 
{ 


system("cls™ ): 


if(wint == 3) 
printf("““RAISED COSINE ARRAY AMPLITUDE WEIGHTS\n\n"): 
else 


printf’ NRAISED COSINE SQUARED ARRAY AMPLITUDE WEIGHTS\o\n"): 
pnntf('NN NUMBER OF ELEMENTS <%d>\n ",m); 
repeats: pnotf("ENTER THE FRACTIONAL PEDESTAL HEIGHT (0.0 - 1.0)\n"): 
pnntf("NOTE A PEDESTAL HEIGHT OF 1.0 RETURNS A RECT WINDOWo"); 
printf“ HEIGHT: " height): 
getsi inline ): 
Stofa(inline. &height,&ncnt.1 ): 
for(j=1.j<=m:j+4+) 


proj = 1.-G-1)*2.(m-1.): 
if(wint == 3) /* raised cos weight calculation */ 
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win{j] = height + (.5-height/2.)*cos(proj*PI); 
else /* raised cos squared weights */ 
Witty} = beigin + (1.-height)*cos(proj*PI/2.)*cos(proj*PI/2.); 


} 
if(2* (m/2) != m) 


num = (m-1.)/2.41.; /* odd number of elements */ 
win[(int) num] = win{(int) num]/2.0; 


} 
printf AWN VERIFY CALCULATED WEIGHT VALUES\\n"); 
for(j=1.j<=m:j++) 


{ 
printf(“’ ELEMENT %d WEIGHT : %8.4f\n "j,win{j]): 


} 
printf("ARE THESE VALUES CORRECT <y/n>? "); 
if(ch = getche() == 'n’) 

{ 


goto repeat3; /* reenter pedestal fractional height */ 


} /* close cosine amplitude weight calculations */ 
if(wint == 5) /* Taylor weight factor calculations */ 
{ 
system("cls"): 
printf(’NTAYLOR FUNCTION ARRAY AMPLITUDE WEIGHTS‘\n"); 
printf(NN’ NUMBER OF ELEMENTS <%d>‘n ",m); 
repeats: printf("ENTER THE TAYLOR PARAMETER (0 - 6)\n"): 
printf "NOTE LARGER VALUES GIVE SMALL SIDELOBES AND WIDE MAIN LOBE\n"); 
printf("(“t TAYLOR PARAMETER: ",height): /* REUSE HEIGHT variable*/ 
gets(inline ): 
Stofa(inline. &height.&ncnt,1): 
for(j=lij<=mij++) 


proj = 1.-(j-1)*2./¢m-1.): 
num = height*sqrt( 1-proj* proj): 
win{j]=bessi0(num): /* weight = mod bessel order 0 */ 


if(2* (m/2) '= m) 
{ 
num = (m-J.)/2.+1.: /* odd number of elements */ 
win[(int) num] = win{(int) num}/2.0: 


} 
printf( \“WNVERIFY CALCULATED WEIGHT VALUES\n\n"); 
for(j=l.j<=m:j++) 


printf“ ELEMENT %d WEIGHT : %8.4fn "j,win{j]); 


| 
printf("ARE THESE VALUES CORRECT <y/n>? "); 
if(ch = getche() = ‘n') 

{ 


goto repeatS; /* reenter Taylor Parameter */ 


} /* close Taylor function weight calculation loop */ 
if(wint == 6) /* Dolph-Chebyshev weight calculations */ 
{ 
system("cls”’): 
printf “DOLPH-CHEB YSHEV FUNCTION ARRAY AMPLITUDE WEIGHTS\o\n"); 
printf NN NUMBER OF ELEMENTS <%d>\n ".m): 
repeat: printf("ENTER THE DESIRED SIDELOBE LEVEL "): 
gets(inline ). 
Stofa(inline.&height.&ncnt,1): /* reuse variable again */ 
/* even or odd # elements */ 
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flag = 0: 
if(2* (m/2) == m) 

flag = 1. /* even number of elements */ 
if(flag == 1) 

num =m/2.; 
else 

num = (m-£.)/2.; 
k=ceil(num); 
a=dvector(0,m):; 
/* Compute chebyshev polynomial */ 
side=pow(L0.,(height/20.)); 
cheb1=((log(side+sqrt(side*side-1)))/(m-1)); 
cheb=.5*(exp(cheb1 )+exp(-cheb1)); 
angle=1.-1./(cheb*cheb): 
for(i=]:i<=k-flag:i++) 

{ 

fact=1.; 

for(j=l.j<=i-1.j++) 


den = (float)j*(m-1-2*i+j)/((i-j)*(i+1-j)): 
fact=fact *angle *den+ 1: 


a[k-iJ=(m-1)*angle* fact: 
} 
a{k]=1L.: 
norm=a|[ flag]: 
if(flag == 0) 
win[k+1]=1: 
for(j=1:j<=k:j++) 


win{j]=a{k+1-j]/norm:; 
win[m+1 -jJ=win{j]; 


if(2* (m/2) != m) 
{ 


num = (m-1.)/2.+1.; /* odd number of elements */ 
win[(int) num) = win[(int) num ]/2.0; 


printf \“WN VERIFY CALCULATED WEIGHT VALUES\o\n"): 
for(j=1:j<=m j++) 


printf‘ ELEMENT %d WEIGHT : %8.4f\n "j,win[j]): 


} 
prnii("ARE THESE VALUES CORRECT <y/n>? "); 
if(ch = getche() == ‘n’) 


goto repeat6; /* reenter sidelobe levels */ 
} 


free_dvector(a,0.num): 
} /* close Chebyshev weight loop */ 


if(wint > 1) /* control loop to calculate beam non rect window */ 
{ 
for(i=0;i<=720;i++) /* rezero ftheta vector */ 

{ 

ftheta[iJ=0.0: 


/* even or odd # elements */ 
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flag = 
if(2* aoe == m) 
{ 


flag = 1; /* even number of elements */ 


} 
num =m/2.; 
} 

else 


num =(m-1.)/2.4+1.; 
} 


norm = 0.: 
for(i=0:i<=360:i++) 
{ 


for(j=1;j<=(int)num;j++) 


if(flag==()) 
fact = (j-1.): 
if( flag==1) 
fact = (j-.5): 
weight = win({(int) num L-j]: 
ftheta[iJ= =ftheta[i]+weight*cos(fact*2.*PI*d/lambda*(cos(i*P1/360.)-cos(theta))); 


ftheta{i]=fabs(ftheta[i]): 
if(ftheta{i] > norm) 
nom = ftheta{[i]: 


} 
for(i=361 .i<=720:14++) 
ftheta[i] = ftheta[720-i]; 
for(i=0:1<=720; i++) 


ftheta[i]=ftheta[i]/norm: 
ftheta[i]=ftheta[i]* ftheta[i]; 
if(ftheta[i}<= .0000001) 
ftheta[i]=.0000001: 
G[i]=10*log10(ftheta[i]}): 


] 

*di=dir( theta); 

*bw=dbdown(G,thetai); 

} /* close beam calc control loop */ 

} /* close beam function */ 

POISE IOIGI III OIG SI IIOIOIOIOIOI IOI IOIS IOI GOI GIOI IU II IOIIOI GI I IOIOIOI IOI IOI I IOI IOI ROI BOtoi a 
* FUNCTION BESSIO(X) CALCULATES MODIFIED BESSEL FUNCTIONS OF ZERO ORDER * 
* ADAPTED FROM NUMERICAL RECIPES IN C BY PRESS ET. AL. 


SHASAA SS EEEESAedeSRRRERNGEHOERMAHOHEHSSHEAETESORMESOESORSORENEEEOEEERH/ 
float bessi0(x) 
float x; 


{ 


float ax, ans: 
double y: 
if((ax = fabs(x)) < 3.75) 


{ 
y=x/3.75; 
yeny: 
ans=1 (+y*(3.5156229+y*(3.0899424+y*(1.2067492+y*(0.2659732+y* 
(0.360768e-1+y*0.4581 3e-2))))): 
} 


else 


{ 
y=3.75/ax: 
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ans=(exp(ax)/sqrt(ax))*(0.39894228+y*(0.1328592e-1+y*(0.2253 19e-2 
+y*(-0.157565e-2+y*(0.91628 Le-2+y*(-0.2057706e-1+y *(0.2635537e-1 
+y*(-0.1647633e-1+y*0.392377e-2)))))))); 
} 
retum ans; 
/* close modified bessel function bessi0() */ 
POGGIO OIG OI GOIOI IOI IOIOIOO IOI I IOI IDI IOIOIOI OIG toi II I IOI soi io ioioi toi tok 
* CALCULATION OF DIRECTIVITY BY NUMERICAL INTEGRATION OF ARRAY 
* FACTOR USING EXTENDED SIMPSON’S RULE 
seesansiatssieetehonastseeneusaeeenesaasonenenensesesensannsseneaneeneesse/ 
float dir(ftheta) 
float ftheta[]: 
{ 


float s; 

int i; 

s=0.0; 
for(i=0;i<=7 1 8;i=1+2) 


{ 
s=8+1.0/(3.0)*.5/360.*(ftheta[i}+4.0* ftheta[i+1]+ftheta[i+2]); 


} 
s=10*logl0(1./s): 
retum s; 


} 


PAF I IOI IOI III IIOIIOIO IOC IOI ISO IOIOIOI SI i iS i loi mai i ioioioi i i daiaieiok 


* FUNCTION TO CALCULATE 3 dB MAIN LOBE BEAMWIDTH * 
BEBO BICC IOI II OI OO TO I IOI I i i ttt 
float dbdown(G.thetai) 

float G[}: 

int thetai: 


float s,sb: 
s=0.0: 
for(i=1:i<=90;i1++) 


sb=s: 
s=G[2*thetai ]-G[(2*thetai)+i]; 
if(fabs(s)>=3.) 

retum (i- 1)+(3-sb)/(s-sb); 


} 
printf(" 3 dB down main lobe width exceeds 90 degrees\n"); 
printf("Press any key to continue"): 
getch(): 
} 
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SO oi ioioiokok ok i lokokokok ik oiioko kkk ek a EEE EEE EEE EE EE 


* PLOT.C A FUNCTION TO PLOT THE CURVES OF NORMALIZED DIRECTIONAL BEAM * 
* PATTERN CALCULATED IN BEAM.C 

Soushuaucaver<escounues 13 se040i exe sa nseoutuauatvasseesvecesvcun euniasueee] 
/* PREPROCESSOR DIRECTIVES */ 

#include <stdio.a> 

#include<graph.h> 

#include <math.h> 

#include <stdlib.h> 

#include <string.h> 

#include “worlddr.h" 

#include "segraph.h" 

/* EXTERNAL VARIABLES DEFINED */ 

float *G; /* pointer to vector of array gain values */ 

int mode: /* integer to select screen mode single vs. dual graph */ 

int num: /* select which array to plot */ 


int m1,m2: /* Number of elements in Array 1, 2 */ 
float f1,f2; /* Signal frequency in Hz areray 1,2 */ 
float di ,d2; /* Interelement spacing in feet array 1,2 */ 
float thetal ,theta2; /* steer "look" angle broadside = 90 */ 
float lambda! .lambda2; /* signal wavelength =) 

float c: /* sound speed in fps */ 

int wint! ,wint2: /* window type */ 

float *G1, *G2; /* pointer for unit offset vector gain */ 
float *wint,*win2; {* pointer for vector of f amplitude tapers */ 
int i. j, k; /* loop counters 

int init: /* flag for initialization of ae windows */ 

float dil ,di2: 


float bw! ,bw?; 
float vmode,pm; 
/* FUNCTION PROTOTYPES */ 
void plot(); 
void plot(mode, num , wint, G) 
float G[]): 
int mode, num, wint; 
{ 
int i, maxX, maxY,plt.data; 
float angle(72!}: — /* vector of x axis data, degrees */ 
char tempstr[20}]: 
cher xstr[20], choice; 
for(i=0:i<=720:1++) 


{ 
angle[i}=i/2.: 
} 


GetMaxCoords(&maxX,&maxY): 

DefGraphWindow(maxX/3, !,maxX,max Y/2,8); 

DefGraph Window(maxX/3,max Y/2,maxX,maxY,10); 
DefGraphWindow(0,0,Round(maxX/3.0),maxY/2.,7); 
DefGraphWindow(0,Round(max Y/2.),Round(maxX/3.),maxY/1.2,9); 
DefGraphWindow(0,Round(max Y/1.2),Round(maxX/3.),maxY,1); 
SetWin2PlotRatio(8,.2,0.14,0.05,0.14): 
SetWin2PlotRatio(10,.2,0.14,0.05,0.14); 
SetWin2PlotRatio(7,0,0,0,0); 

SetWin2PlotRatio(9,0,0,0.0): 

SetWin2PlotRatio(1,0,0,0,0): 


if(num == 1) 
{ 
SetCurrent Window(8): 


ClearWindow(); 
SetAxesType(0,0), 
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if(vmode!=12.) 

SetPlotBackground(3); 
SelectColor( 15); 
AutoAxes(angle,G,720,1); 
TitleWindow("ARRAY #1 BEAM PATTERN"), 
TitleXAxis("Theta (degrees)"): 
Title YAxis("G(Theta) (db)"); 
BorderCurrent Window(4); 
SetCurrent Window(7); 
ClearWindow(); 
if(vmode!=12.) 

SetPlotBackground(3); 
BorderCurrent Window(4); 
SelectColor(15); 
setlinestyleXX(0,0,1); 
TitleWindow("ARRAY #1"); 
settextjustifyXX(0,0); 
sprintf(xstr,"%6.0f" £1); 
strcpy(tempstr,"Freq (Hz): "); 
strcat(tempstr, xstr ); 
LabelGraph Window(20,880,tempstr,0,0); 
sprintf(xstr."%6.0f",c); 
strcpy(tempstr,"Celerity (fps): ”); 
strcat(tempstr, xstr ): 
LabelGraph Window(20,830,tempstr,0,0); 
sprintf(xstr,"%3d".m1): 
strcpy(tempstr,"# elements: "): 
strcat(tempstr, xstr ): 
Label Graph Window(20,780,tempstr,0,0): 
sprintf(xstr,"%5.2f" 1): 
strcpy(tempstr, "spacing (ft): "); 
strcat(tempstr, xstr ); 
LabelGraph Window(20,730,tempstr,0,0); 
sprintf(xstr,"%4.2f" lambda} ); 
strcpy(tempstr," Wavelength (ft): "); 
strcat(tempstr, xstr ); 
LabelGraph Window(20,680,tempstr,0,0); 
sprintf(xstr,"%4. 1f" theta 1); 
strcpy(tempstr,”"Look Angle (deg): "); 
strcat(tempstr, xstr ); 
LabelGraph Window(20,630,tempstr,0,0); 
switch(wint) 


case |: 
strcpy(tempstr,"Window: RECTANGULAR"); 
LabelGraphWindow(20,580,tempstr,0,0); 
break; 

case 2: 
strcpy(tempstr,"Window: USER"); 
LabelGraphWindow(20,580,tempstr,0,0); 
break: 

case 3: 
strcpy(tempstr,"Window: COSINE"); 
LabelGraphWindow(20,580,tempstr,0,0): 
break; 

case 4: 
strcpy(tempstr," window: COSINE SQUARED"); 
LabelGraphWindow(20,580,tempstr,0,0); 
break; 

case 5: 
strcpy(tempstr,”"Window: TAYLOR"); 
LabelGraphWindow(20,580,tempstr,0,0); 
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break; 

case 6: 
strcpy(tempstr,"Window: CHEB YSHEV"); 
LabelGraphWindow(20,580,tempstr,0,0); 
break; 


} 
sprintf(xstr,"%5.3f" dil): 
strcpy(tempstr,"Directivity(dB): "); 
Strcat(tempstr, xstr ); 
LabelGraphWindow(20,530,tempstr,0,0); 
sprintf(xstr,"%2. If" ,bw 1); 
strcpy(tempstr,"3 dB beam width: "); 
strcat(tempstr, xstr ); 
LabelGraph Window(20,480,tempstr,0,0): 


SetCurrent Window(8); 
ClearGraphi): 
if(vmode!=12.) 
SetPlotBack ground(3); 
SelectColor(4): 
DrawGnid(6): 
LinePlotData(angle,G,721,15,0); /* command to plot array #1 G(theta) */ 
} /* close num = | */ 
else 


{ 

SetCurrent Window(10); 
ClearWindow(): 
SetAxesType(0,0); 
if(vmode!=12.) 

SetPlotBack ground( 14); 
SelectColor(15); 
AutoAxes(angle,G,720,1); 

Title Window("ARRAY #2 BEAM PATTERN’); 
TitleXAxis("Theta (degrees)"); 
Title Y Axis("G(Theta) (db)"); 

BorderCurrent Window(3):; 
SetCurrentWindow(9); 

ClearWindow(); 

if(vmode!=12.) 

SetPlotBack ground(14): 

BorderCurrent Window(3): 
SelectColor(1): 

setlinestyleXX(0,0,1): 

Title Window("ARRAY #2"); 
SettextjustifyXX(0,0); 
sprintf(xstr,"%6.0f" £2); 
strcpy(tempstr," Freq (Hz): "); 
Strcat(tempstr, xstr ); 

LabelGraph Window(20,850,tempstr,0,0); 
sprintf(xstr,"%6.0f',c); 
strcpy(tempstr,"Celerity (fps): "); 
strcat(tempstr, xstr ): 

LabelGraph Window(20,780,tempstr,0,0); 
sprintf(xstr,"%3d",m2); 
strcpy(tempstr,"# elements: "); 
strcat(tempsh, xstr ): 
LabelGraph Window (20,710.tempstr,0,0); 
sprintf(xstr,"%5.2f",d2); 
strcpy(tempstr,"spacing (ft): "): 
strcat(tempstr, xstr ); 

LabelGraph Window(20,640,tempstr,0,0): 
sprint{(xstr,"%4.2f" lambda2); 
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strcpy(tempstr," Wavelength (ft): "); 
Strcat(tempstr, xstr ); 
LabelGraphWindow(20,570,tempstr,0,0); 
sprintf(xstr,"%4. 1f" ,theta2); 
strcpy(tempstr,"Look Angle (deg): "); 
strcat(tempstr, xstr ); 
LabelGraph Window (20,490,tempstr,0,0); 
switch(wint) 
{ 
case |: 
strcpy(tempstr,"Window: RECTANGULAR"); 
LabelGraphWindow(20,4 10,tempstr,0,0); 
break: 
case 2: 
strcpy(tempstr,"Window: USER"); 
LabelGraphWindow(20,410,tempstr,0,0); 
break: 
case 3: 
strcpy(tempstr,"Window: COSINE"); 
LabelGraphWindow(20,410,tempstr,0,0); 
break: 
case 4: 
strcpy(tempstr,” Window: COSINE SQUARED"): 
LabelGraphWindow(20,410,tempstr,0,0); 
break: 
case 5: 
strcpy(tempstr,"Window: TAYLOR"): 


LabelGraphW indow(20.410,tempstr,0,0); 
break: 


case 6: 
strcpy(tempstr," Window: CHEB YSHEV"); 
LabelGraphWiadow(20,410,tempstr,0,0); 
break; 
} 
sprintf(xstr,"%5.3f" ,di2); 
strcpy(tempstr, Directivity(dB): "); 
strcat(tempstr, xstr ): 
LabelGraph Window(20,330,tempstr,0,0); 
sprintf(xstr,"%2.1f"' ,bw2): 
strcpy(tempstr,"3 dB beam width: "); 
strcat(tempstr, xstr ): 
LabelGraphWindow(20.250,tempstr,0,0): 
SetCurrent Window(10): 
ClearGraph(): 
if(vmode!=12.) 
SetPlotBackground(14); 
SelectColor(2); 
DrawGrid(6); 
LinePlotData(angle,G,720,1,0); /* COMMAND TO PLOT ARRAY 2 DATA */ 
} /* else */ 
} /* close plot function */ 











POR BIOI IOIOI IOI I IO IOI II IIS ISIC OI IO IOI IDI pi toloi i tok i toi i tok 


* —QUTPUT.C CONTROLS ARRAY DATA OUTPUT TO SELECTED FILE 
* OR TO PRINTER 

alaaaadinta didi dlaadaelaaeialatelae dedicated clea diaaeadialaceclac acacia tila ceded dla dndhadhcicethatdilaatadhaca eth iceddtd | 
/* PREPROCESSOR DIRECTIVES */ 

/* PREPROCESSOR DIRECTIVES */ 

#include <stdio.h> 

#include <setjmp.h> 

#include <math.h> 

#include <stdlib.h> 

#include <graph.h> 

#include <string.h> 

#include "worlddr.h" 

#include "se graph.h" 


/* EXTERNAL VARIABLES DEFINED */ 


int m: /* Number of elements in Array 1, 2 */ 
float f: /* Signal frequency in Hz areray 1,2 */ 
float d: /* Interelement spacing in feet array 1,2 */ 
float theta: /* steer "look" angle broadside = 90) */ 
float lambda; /* signal wavelength */ 

float c: /* sound speed in fps */ 

int wint: /* window type */ 

float *G; /* pointer for unit offset vector gain */ 
float *win: /* pointer for vector of amplitude tapers */ 
int i, j. k: /* loop counters od 

int wir: /* amplitude taper type marker */ 
char inline[8 1]; /* data entry */ 

char nent: 

float di: 

float bw: 

void output(): 

float pm: 


void output(f,c,m,dlambda, theta, wint,di,ow,win,G) 
float f,c,d.lambda,theta,di,bw; 
float G[}. win[]: 
int m,wint: 
{ 
int error. 
int i.maxX.maxY; 
float angle[721}:  /* vector of x axis data, degrees */ 
char tempstr[ 20]: 
char xstr[20], choice; 
FILE *out: 
for(i=0:1<=720:1++4) 


{ 
angle[i}=i/2.: 
} 


out= fopen(“array.dat”,"a+"); /* open array.dat for appending */ 


| haslacladlachiciaahadiaaiaalaabasiacladeleatealealteciadladlacladhadladladladhaihachaciadiadintaiciaiadasladiadladiadhasbaalachachalashashecheleachedhachshshachahshdhaaahaaashaheeaaedeaeaiel 


* OUTPUT DATA TO FILE ARRAY.DAT CREATED/OPENED ON DEFAULT DRIVE 
BEER BIE OIG GIGS OIC IE CI ISI IG I III i A IOI I I i RR A a i Rota teat ie te / 

fprintf(out, WWWARRAY OUTPUT DATA FILE\o\n"); 

fprintf(out,” Array Design Characteristics:\n"); 

fprintf(out,"“tFrequency NN %6.0f Hz\n",f); 

rintf(out,'Ntcelerity NN%4.0f fps\n",c); 

Ponntf out‘ elementsNN%d\n ",m); 

fprintf(out, NSpacingNN%4.1f ft\n'".d): 

fprintf(out, NwavelengthNN%3.1f ft\n" lambda); 

switch(wint) 
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{ 

case |: 
fprintf(out,, “Window NNRECTn"); 
break; 

case 2: 
fprintf(out, “Window NNNUSER\n"); 
break: 

case 3: 
fprintf(out,, “Window NNNCOSINE\n"); 
break: 

case 4: 
fprintf(out, “Window NNNCOSINE SQUARED\n"): 
break: 

case 5: 
fprintf(out,, “Window: NNNTA YLOR\n"); 
break; 


case 6: 
fprintf(out,"\“Window:NNNCHEB YSHEV\n"); 
break: 


} 
fprintf(out,’\Directivity NN%5.3f dB\n",d): 
fprintf(out.'\wBeam width NN%4.2f degrees\n\n" bw); 
if(wint != 1) 


{ 
fprintf(out," ARRAY WINDOW WEIGHTS\n\n"): 
for(i=l:i<=m:i++) 

fprintf(out,’wElement[%d] Amplitude Weight: %f\n",i,win{[i]); 


} 
fprintf(out," WNORMALIZED DIRECTIONAL GAIN PATTERN VALUES IN dB\n\n"); 
fort i=1:i<=360;i=i+4) 
fprintflout."%d deg: %4.2N%d deg: 4.2NKd deg: %4.2NI%d deg: %4.2fn", 
i,G[2*i].i+1,G[2*(i+1)),i+2,G[2 *(i+2)].i+3,G[2*(i+3))): 
fprintf(out."\f"): 
fclose(out): : 
FhaciathachastaalachaciadiadladiaclasIaalealealacladiacladladieshectaatadiaciadladiadladhadindindinsinalaciadiasiaaladiadiachadlasteelesleslaalaalediediaslasladleslasindlastiladlacieshaslaslasiadlachealeaadah dasa 


PRODUCE PRINTER OUTPUT OF NORMALIZED DIRECTIONAL GAIN * 
RRR REE EEE RE ER EERE REAR EEE EE EEE REE E EEE EEE EE EEE EEER EEE EEE EE EH / 
GetMaxCoords(&maxX,&maxY): 
DefGraphWindow(0,0,Round(maxX/3.),maxY,5); 
DefGraphWindow(Round(maxX/3),0,maxX,maxY,6); 

Set Win2PlotRatio(5,.2,0.14,0.05,0.14): 

SetWin2PlotRatio(6,.2,0.14 9.05,0.14); 

SetCurrent Window(6); 

ClearWindow(); 

SetAxesType(0,0); 

SelectColor( 15): 

AutoAxes(angle,G,720,1); 

TitleWindow("ARRAY BEAM PATTERN"); 

TitleXAxis("Theta (degrees)"); 

Title Y Axis("G(Theta) (db)"): 

BorderCurrent Window(4): 

SelectColor(4): 

DrawGrid(6): 

LinePlotData(angle,G,721,15,0); /* command to plot array #1 G(theta) */ 
SetCurrent Window(5): 

ClearWindow(): 

BorderCurrent Window(4): 

SelectColor( 15): 

setlinestyleXX(0,0,!); 

TitleWindow("ARRAY CHARACTERISTICS"); 

settextjustifyXX(0,0); 

sprintf(xstr,"%6.0f",f); 
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strcpy(tempstr,"Freq (Hz): ")s 
strcat(tempstr, xstr ); 
LabelGraph Window (20,880,tempstr,0,0); 
sprintf(xstr,"%6.0f",c); 
strcpy(tempstr,"Celerity (fps): "); 
strcat(tempst, xstr ); 
LabelGraph Window(20.830,tempstr,0,0); 
sprintt(xstr,"%3d",m); 
strcpy(tempstr, "# elements: "); 
Strcat(tempstr, xstr ); 
LabelGraph Window(20,780,tempstr,0,0); 
sprintf(xstr,"%5.2f'.d); 
strcpy(tempstr, "spacing (ft): "); 
strcat(tempstr, xstr ); 
LabelGraphWindow(20,730,tempstr,0,0); 
sprintf(xstr,"%4.2f" lambda); 
strcpy(tempstr,"Wavelength (ft): "); 
strcat(tempstr, xstr ); 
LabelGraphWindow(20,680,tempstr,0,0); 
sprintf(xstr."%4. If" theta); 
strcpy(tempstr,""Look Angle (deg): "); 
strcat(tempstr, xstr ); 
LabelGraphWindow(20,630,tempstr,0,0): 
switch(wint) 
{ 
case |: 
strcpy(tempstr,"Window: RECTANGULAR’), 
LabelGraphWindo w:20,580,tempstr,0,0); 
break: 
case 2: 
strcpy(tempstr,"Window: USER"); 
LabelGrapi:Window(20,580,tempstr,0,0); 
break: 
case 3: 
strepy(tempstr,” Window: COSINE"), 
LabelGraphWindow(20,580,tempstr,0,0); 
break: 
case 4: 
strcpy(tempstr," Window: COSINE SQUARED"); 
LabelGraphWindow(20,580,tempstr,0,U): 
break: 
case 5: 
strcpy(tempstr."Window: TAYLOR”): 
LabelGraphWindow(20,580,tempstr.0,0): 
break: 
case 6: 
strcpy(tempstr,"Window: CHEB YSHEV"): 
LabelGraphWindow(20.580,tempstr,0,0): 
break: 
} 
sprintf(xstr,"%5.3F',di): 
strcpy(tempstr,"Directivity(dB): "); 
strcat(tempst,, xstr ): 
LabelGraph Window(20,530,tempstr,0,0); 
sprintf(xstr,"%2.11" ,bw): 
strcpy(tempstr,"3 dB beam width: "); 
strcai(tempstr, xstr ): 
LabelGraphWindow(20,480,tempstr,0,0); 
sircpy(tempstr,"ENSURE PRINTER IS READY"), 
LabelGraph Window(20,420,tempstr,0,0): 
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LabelGraphWindow(20,300,"PRESS ANY KEY TO PRINT",0,0); 
choice = getch(); 
ScreenDump((int) pm,0,5,1.1,1.0,0,1,0,&error); 


} 





95 








Appendix B ARRAY Output During Synthesis 
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Array Characteristics 


ARRAY CHARACTERISTICS | ARRAY BEAM PATTERN 


Free (itz>: 2e0e] 
Celerity (fps): aceu 
Se elementa: 8 
spacing CF ED: 1.22 
Wavetength ¢€ft): 2.43 
Look Aanwle Cde@~>: 938.9 
Hindow: RECTANGULAR 
DireuntivityCaB): 

3 dB beam width: 

ENSURE PRINTER 1S READY | 


G(Theta? (db? 


PRESS ANY KEY TO PRINT f 
i 


160 zee 


Theta Cdegrees) 





Frequency (Hz): 2000 
Sound Speed( fps): 4860 
# Elements: 8 
Element Spacing(ft): 1.2 
Wavelength (ft): 2.4 
Directivity (dB): 1.215 
Beamwidth (degrees): 12.77 
Degrees G[9] Degrees G[8] Degrees G[@] Degrees G6) 
1 -70 2 -60.38 3 -53.34 4 -48.35 
5 -44.47 6 -41.31 7 -38.64 8 -36.33 
9 -34,3 10 -32.5 i} -30.87 12 -29.39 
13 -28.05 14 -26.81 15 -25.68 16 -24.64 
17 -23.69 18 -22.81 19 -22.01 20 -21.28 
21 -20.62 22 -20.02 23 -19.5 24 -19.04 
25 - 18.66 26 ~18.35 27 -18.11 28 -17.96 
29 -17.89 30 -17.92 31 -18.06 32 -18.32 
33 -18.72 34 -19,28 35 -20.05 36 -21.07 
37 -22.44 38 -24.32 39 -27.03 40 -31.43 
41 -41.96 42 -38.64 43 -29.95 44 -25.69 
45 -22.9 46 -20.89 47 -19.38 48 -18.24 
49 -17.41 50 -16.83 5] -16.51 52 -16.43 
53 -16.61 54 -17.08 55 -17.9 56 -19.17 
57 -21.1 58 -24.14 59 -29.78 60 -70 
61 -29.28 62 -23.15 63 -19.62 64 -17.23 
65 -15.51 66 -14.28 67 -13.44 68 -12.95 
69 -12.8 70 -13.01 7 -13.65 72 -14.85 
73 - 16.88 74 -20.51 75 -29.13 76 -29,33 
77 -19.03 78 -14.13 79 -10.85 80 -8.41 
81 -6.48 82 -4.92 83 -3.66 84 -2.62 
85 -1.78 86 -1.12 87 -0.63 88 -0.28 
89 -0.07 90 0 91 -0.07 92 -0.28 
93 -0.63 94 -1.12 95 -1.78 96 -2.62 
97 -3.66 98 -4.92 99 -6.48 100 -8.41 
10! -10.85 102 -14.13 103 -19.03 104 -29.33 
105 -29.13 106 -20.51 107 -16.88 108 -14.85 
109 -13.65 110 -13.01 M11 -12.8 112 -12.95 
113 -13.44 114 -14.28 115 -15.51 116 -17.23 
117 -19.62 118 -23.15 119 -29.28 120 -70 
121 -29.78 122 -24.14 123 -21.1 124 -19.17 
125 -17.9 126 -17.08 127 -16.61 128 -16.43 
129 -16.51 130 -16.83 13) -17.41 132 -18.24 
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133 
137 
141 
145 
149 
153 
157 
161 
165 
169 
173 
177 


-19.38 
-29.95 
-27.03 
-20.05 
-18.06 
-18.11 
-19.5 
-22.01 
-25.68 
-30.87 
-38.64 
-53.34 











Array Characteristics 


ARRAY CHARACTERISTICS 


Freq (Hz): 

Celerityu ¢fyps): 

WH elements: 2 
spacing Ct>?: 1.22 
Wavelength (ft): 8.97 
Look Angie €deg>: 383.03 
Window: RECTANGULAR 
DirectivitycdB)d: 8.5387 
A AR heam width: SS. 1 


ENSURE PRINTER 15 READY 


G(Thetad (db) 


PRESS ANY KEY TO PRINT 


Frequency (Hz): 5000 
Sound Speedifps): 4860 
# Elements: 8 
Eiement Spacing(ft): 1.2 
Wavelength (ft): 1 
Directivity (dB): 1.215 


Beamwidth 5.09 
(degrees): 

Degrees G[9] Degrees G{e] 
1 -61.45 2 -49,39 
5 -33.39 6 -30.19 
9 -23.08 10 -21.26 
13 -16.92 14 -15.81 
17 -13.44 18 -13.02 
21 -13.27 22 -14.12 
25 -24.34 26 -38.05 
29 -10.18 30 -7.39 
33 -2.22 34 -4.22 
37 0 38 -0.2 
41 3.11 42 -5.17 
45 -22.21 46 -25.69 
49 -12.8 50 -13.54 
53 -39.87 54 -23.98 
57 -16.7 58 -18.74 
61 -24.33 62 -19.47 
65 -21.74 66 -31.61 
69 -18.25 70 -18.11 
73 -29.75 74 -20.17 
77 -18.85 78 -27.49 
81 -13.45 82 -12.89 
85 -16.7 86 -8.5 
89 -0.43 90 0 
93 -4,28 94 -8.5 
97 -15.28 98 -12.89 

101 -25.64 102 -27.49 
105 -16.99 106 -20.17 
109 -20.54 110 “18.11 
113 -28.9 114 -31.61 
117 -17.93 118 -19.47 
121 -23.98 122 -18.74 
125 -18.5 126 -23.98 


ARRAY BEAM PATTERN 


i 


Bieatantwiert te 
TNC 
3 ML oI 
Boevesiiics 


1898 2a8 


Theta (degrees) 





Degrees G[@] Degrees G8} 
3 


-42.33 4 -37.3 
7 -27.48 8 -25.14 
il -19.64 12 -18.2 
15 -14.85 16 -14.06 
19 -12.81 20 -12.87 
23 -15.67 24 -18.44 
27 -20.03 28 -14.01 
31 -5.23 32 -3.54 
35 -0.52 36 -0.11 
39 -0.75 40 -1.69 
43 -8.15 44 -12.8 
47 -16.49 48 -13.57 
51 -15.9 52 -21 
55 -18.5 56 -16.61 
59 -23.98 60 -70 
63 -17.93 64 -18.56 
67 -28.9 68 -20.92 
71 -20.54 72 -28.57 
75 -16.99 76 -16.55 
79 -25.64 80 -16.71 
83 -15.28 84 -27.05 
87 -4.28 88 -1.79 
91 -0.43 92 -1.79 
95 -16.7 96 -27.05 
99 -13.45 100 -16.71 
103 -18.85 104 -16.55 
107 -29.75 108 -28.57 
fl -18.25 112 -20.92 
115 -21.74 116 -18.56 
119 -24,33 120 -70 
123 -16.7 124 -16.61 
127 -39.87 128 -21 
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129 
133 
137 
141 
145 
149 
183 
157 
161 
165 
169 
173 
177 


-15.9 
-16.49 
-8.15 
-0.75 
-0.52 
-5.23 
-20.03 
-15.67 
-12.81 
-14.85 
-19.64 
-27.48 
-42.33 


130 
134 
138 
142 
146 
150 
154 
158 
162 
166 
170 
174 
178 


-13.54 
-25.69 
-5.17 
-0.2 
-1.22 
-7.39 
-38.05 
-14.12 
-13.02 
-15.81 
-21.26 
-30.19 
-49.39 


131 
135 
139 
143 
147 
151 
155 
159 
163 
167 
171 
175 
179 


100 


-13.57 
-12.8 
-1.69 
-0.11 
-3.54 

-14.01 

-18.44 

-12.87 

-14.06 
-18.2 

-25.14 
-37.3 
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Array Characteristics 






ARRAY BEAM PATTERN 


ARRAY CHARACTERISTICS 
Freq (Hz>?: 29690 
Celerityw (fps): 4860 
Ww elements: is 


wrpacing CFt?: @.36 
Wavelength Cf: 2.43 
Look Angle Cdeg>: 90.8 
Hindow: RECTANGULAR 
DirectivitycdB>?: 18.439 
3 A&B heam width: 14.2 
ENSURE PRINTER 18 READY 


G(Theta) (dh) 


PRESS ANY KEY TO PRINT 


| Theta (degrees) 








Frequency (Hz): 2000 
Sound Speed(fps): 4860 
# Elements: 15 
Element Spacing(ft): 0.6 
Wavelength (ft): 2.4 
Directivity (dB): 0.581 
Beamwidth (degrees): 14.2 
Degrees G[8] Degrees G[8] Degrees G[8] Degrees G8] 
1 -20.49 2 -20.48 3 -20.45 4 -20.42 
5 -20.38 6 -20.33 7 -20.28 8 -20.23 
9 -20.18 10 -20.13 il -20.08 12 -20.04 
13 -20.02 14 -20.01 15 -20.02 16 -20.05 
17 -20.11 18 -20.21 19 -20.35 20 -20.53 
21 -20.78 22 -21.09 23 -21.48 24 -21.97 
25 -22.57 26 -23.32 27 -24.26 28 -25.46 
29 -27 30 -29.09 31 -32.12 32 -37.3 
33 -54.6 34 -39.57 35 -32.64 36 -28.76 
37 -26.07 38 -24.03 39 -22.41 40 -21.11 
4\ -20.04 42 -19.19 43 -18.51 44 -18 
45 -17.66 46 -17.47 47 -17.46 48 -17.63 
49 -18 50 -18.63 51 -19.56 52 -20.89 
53 -22.84 54 -25.82 55 -31.17 56 -54.24 
57 -31.91 58 -25.39 59 -21.67 60 -19.11 
61 -17.22 62 -15.79 63 -14.71 64 -13.93 
65 -13.41 66 -13.16 67 -13.18 68 -13.5 
69 -14.18 70 -15.35 71 -17.22 72 -20.36 
73 -26.8 74 -38.03 75 -22.15 76 -16.48 
77 -12.88 78 -10.23 79 -8.15 80 -6.45 
81 -5.05 82 -3.88 83 -2.91 84 -2.1 
85 -1.44 86 -0.91 87 -0.51 88 -0.22 
89 -0.06 90 0 91 -0.06 92 -0.22 
93 -0.51 94 -0.91 95 -1.44 96 -2.1 
97 -2.91 98 -3.88 99 -5.05 100 -6.45 
101 -8.15 102 -10.23 103 -12.88 104 -16.48 
105 -22.15 106 -38.03 107 -26.8 108 -20.36 
109 -17.22 110 -15.35 iil -14.18 112 -13.5 
113 -13.18 114 -13.16 115 -13.41 116 -13.93 
17 -14.71 118 -15.79 119 -17.22 120 -19.11 
121 -21.67 122 -25.39 123 aay at 124 -54.24 
12 -3017 126 -25.82 127 -22.84 128 -20.89 
129 -19.56 130 -18.63 131 -18 132 -17.63 











133 
137 
14] 
145 
149 
153 
457 
161 
165 
169 
173 
177 





-17.46 
-18.51 
-22.41 
-32.64 
-32.12 
-24.26 
-21.48 
-20.35 
-20.02 
-20.08 
-20.28 
-20.45 


134 
138 
142 
146 
150 
154 
158 
162 
166 
170 
174 
178 


-17.47 
-19.19 
-24.03 
-39.57 
-29.09 
-23.32 
-21.09 
-20.21 
-20.01 
-20.13 
-20.33 
-20.48 


135 
139 
143 
147 
151 
155 
159 
163 
167 
171 
175 
179 
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-17.66 
-20.04 
-26.07 
-54.6 
-27 
-22.57 
-20.78 
-20.11 
-20.02 
-20.18 
-20.38 
-20.49 















Freq (tx): 





Celerity (fps): 
Mw element«: 


spaoing (ft): 





Wavelength Cftd: 





Hindow: 
DirectivitywC dB): 
3 aB heam width: 


ENSURE PRINTER (3 READY 


PRESS ANY KEY TO PRINT 


ARRAY CHARACTERISTICS 


Look Angle Cdeaq>: 
RECTANCULAR 
14.367 








Array Characteristics 


C(Theta) (db?) 





Frequency (Hz): 5000 
Sound Speed( fps): 4860 
# Elements: 15 
Element Spacing(ft): : 6 
Wavelength (ft): 

Directivity (dB): 0. 581 


Beamwidth (degrees): 5.64 


Degrees ae Degrees 
I 2 


-40. 

S -35. 68 6 

9 -29.93 10 
13 ~-25.7 14 
17 -23.5 18 
21 -24.21 22 
25 -32.2 26 
29 -29.61 30 
33 -23.53 34 
37 -30.12 38 
Al -26.74 42 
45 -24.49 46 
49 -31.5 50 
53 -23.16 54 
57 -30.76 58 
61 -22.57 62 
65 -24.69 66 
69 -23.65 70 
73 -18.49 74 
77 -42.2 78 
81 -13.17 82 
85 -11.69 86 
89 -0.35 90 
93 -3.39 94 
97 -21.56 98 
101 -16.44 102 
105 -18.4 106 
109 -30.42 110 
413 -20.03 114 
117 -34.65 118 
121 -22.4 122 
125 -29.93 126 
129 -23.68 130 





G[9] 
-39.74 
-34.16 

-28.7 
-24.94 
-23.34 
-25.12 
-39.29 
-26.75 
-23.75 
-38.14 
-24.49 
-27.02 
-26.17 
-25.14 
-24.92 
-25.67 
-21.23 
-33.07 
-17.44 
-22.69 
-14.86 

-6.52 

0 

-6.52 
-14.86 
-22.69 
-17.44 
-33.07 
-21.23 
-25.67 
-24,92 
-25.14 
-26.17 






if il 


| 


ine= 


Degrees 
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RRRAY BEAM PRITERN 


19e¢ 


G[8] 
-38.58 
-32.67 
-27.58 
-24.31 
-23.38 
-26.53 

-48.4 
-24.96 
-24.69 
-41.93 
-23.46 
-32.58 
-23.68 
-29.93 

-22.4 
-34.65 
-20.03 
-30.42 

-18.4 
-16.44 
-21.56 

-3.39 

-0.35 
-11.69 
-13.17 
-42.2 
-18.49 
-23.65 
-24.69 
-22.57 
-30.76 
-23.16 

-31.5 


Theta (degrees> 


Degrees 









G8] 
-37.18 
-31.25 
-26.58 
-23.83 
-23.65 
-28.69 
-34.55 
-23.93 
-26.59 
-31.01 
-23.44 
-58.01 
-22.76 
-52.84 
-21.67 

-33.7 
-20.65 
-21.83 
-22.36 
-13.76 
-23.36 

-1.44 

-1.44 
-23.36 
-13.76 
-22.36 
-21.83 
-20.65 

-33.7 
-21.67 
-52.84 
-22.76 
-58.01 





133 
137 
141 
145 
149 
153 
157 
161 
165 
169 
173 
177 


-32.58 
-23.46 
-41.93 
-24.69 
-24.96 
-48.4 

-26.53 
-23.38 
-24.31 
-27.58 
-32.67 
-38.58 


134 
138 
142 
146 
150 
154 
158 
162 
166 
170 
174 
178 


-27.02 
-24.49 
-38.14 
-23.75 
-26.75 
-39.29 
-25.12 
-23.34 
-24.94 
-28.7 
-34.16 
-39.74 


135 
139 
143 
147 
151 
155 
159 
163 
167 
171 
175 
179 


104 


-24.49 
-26.74 
-30.12 
-23.53 
-29.61 
-32.2 
-24.21 
-23.5 
-25.7 
-29.93 
-35.68 
-40.52 





-23.44 
31.01 
-26.59 
-23.93 
-34.55 
-28.69 
-23.65 
-23.83 
-26.58 
-31.25 
-37.18 
-40.8 








Array Characteristics 


ARRAY CHARACTERISTICS Q@RRAY BEAM PATIEZERN 


Freq (siz? 
Celerity Cfps): 
4 element«: is 


spaacing (ft): 9.38 





Havelength €ft>: 2.43 
Look Anqle (deq>: 68.9 
Hindow: RECTANGULAR 


FECL 


I Le te 





so epee 





wi 
eee aia 
ibs 1m j 
Directivitycap>: 9.023 | 2 COC ET 
3 AB heam width: 15.8 : i Hy 
ENSURE PRINTEK IS READY « La 
PRESS ANY KEY TO PRINT ‘ — 
|_| AEC 
= a 
Theta (Cdegrees> 
Frequency (Hz): 2000 
Sound Speed(fps): 4860 
# Elements: 15 
Element Spacing(ft): 0.6 
Wavelength (ft): 2.4 
Directivity (dB): 0.581 


Beamwidth (degrees): 15.78 
Degrees G[6] Degrees G[8] Degrees G[8] Degrees G@] 
| -19.09 2 -19.03 3 -18.92 4 


-18.77 

5 -18.59 6 -18.37 7 -18.12 8 -17.84 
9 -17.54 10 -17.22 11 -16.89 12 -16.55 
13 -16.2 14 -15.86 1S “15.51 16 -15.18 
17 -14.85 18 -14.54 19 -14.25 20 -13.98 
21 -13.74 22 -13.53 23 -13.36 24 -13.23 
25 -13.15 26 -13.13 27 -13.17 28 -13.29 
29 -13.5 30 -13.81 31 -14.25 32 -14.85 
33 -15.66 34 -16.75 35 -18.23 36 -20.36 
37 -23.67 38 -30.11 39 -43.11 40 -26.18 
41 -20.44 42 -16.8 43 -14.09 44 -11.92 
45 -10.12 46 -8.58 47 -7.24 48 -6.06 
49 -5.03 50 -4.11 51 -3.31 52 -2.6 
53 -1.98 54 -1.46 55 -1.01 56 -0.65 
57 -0.37 58 -0.16 59 -0.04 60 0 
61 -0.04 62 -0.17 63 -0.39 64 -0.71 
65 -1.12 66 -1.65 67 -2.3 68 -3.09 
69 -4.04 70 -5.17 71 -6.52 72 -8.17 
73 -10.19 74 -12.77 75 -16.29 76 -21.82 
77 -36.64 78 -26.98 79 -20.32 80 -17.11 
81 -15.21 82 -14.05 83 -13.4 84 -13.14 
85 -13.22 86 -13.63 87 -14.36 88 -15.45 
89 -16.99 90 -19.11 91 -22.17 92 -27.11 
93 -39.33 94 -33.57 95 -25.99 96 -22.4 
97 -20.23 98 -18.84 99 -17.99 100 -17.55 
101 -17.45 102 -17.67 103 -18.22 104 -19.11 
105 -20.4 106 -22.22 107 -24.8 108 -28.76 
109 -36.48 110 -44.28 111 -31.64 112 -27 
113 -24.,3 114 -22.53 115 -21.35 116 -20.58 
117 -20.15 118 -20.01 119 -20.13 120 -20.5 
121 -21.13 122 -22.06 123 -23.33 124 -25.04 
125 -27.38 126 -30.8 127 -36.66 128 ~-64.1 


129 -37.72 130 -31.73 131 -28.44 132 ~-26.27 








133 
137 
141 
145 
149 
153 
157 
161 
165 
169 
173 
177 











-24.74 
-21.9 
-22.03 
-24.39 
-29.73 
-48.48 
-33.44 
-27.88 
-25.45 
-24.24 
-23.63 
-23.37 


134 
138 
142 
146 
150 
154 
158 
162 
166 
170 
174 
178 


~23.62 
-21.71 
-22.41 
-25.37 
-31.94 
-50.51 
~31.49 
-27.1 
~25.07 
~24.04 
-23.54 
+23.34 


135 
139 
143 
147 
151 
155 
159 
163 
167 
171 
175 
179 
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-22.82 
-21.68 
-22.93 
~-26.55 
-34.91 
-40.51 
-30 
-26.45 
-24.74 
-23.88 
-23.47 
-23.32 


-22.26 
-21.79 
-23.58 
-27.98 
-39.36 
-36.17 
-28.83 
-25.91 
-24.47 
-23.75 
-23.41 
-23.32 











| 
t 


ARRAY CHARACTERISTICS 


Freres ¢Hz>: 2eee 
Celerityu (fps: 4860 
a elements: as 
sapaeucing CHD: @.36 
Wavelength (ft): 2.43 


Look Angle Cdew)d: 982.0 
wHieduw: CHEBYSHEY 
DirectivitucadaB): 9.828 
3 @€B heam width: 147.9 


ENSURE PRINTER 18 READY 


PRESS ANY KEY TO PRINT 


Array Characteristics 


a 
a 
=. 
~ 
A“ 

« 
- 

» 
« 
- 
v 


ARRAY BEAM PATTERN 


Theta (degrees) 





Frequency (Hz): 2000 
Sound Speed( fps): 4860 
# Elements: 15 
Element Spacing(ft): 0.6 
Wavelength (ft): 2.4 
Directivity (dB): 
Beamwidth (degrees): 0.581 
17.89 
Element Amplitude Weight 
1 0.281302 
2 0.326266 
3 0.474929 
4 0.630153 
5 0.775678 
6 0.894747 
7 0.972813 
8 0.5 
9 0.972813 
10 0.894747 
im 0.775678 
12 0.630153 
13 0.474929 
14 0.326266 
15 0.281302 
Degrees G(é@] Degrees 
1 -30.71 2 
5 -30.56 4 
9 -30.29 10 
13 -30.05 14 
17 -30.05 18 
21 -30.62 22 
25 -32.29 26 
29 -36.39 30 
33 -53.23 34 
37 -37.86 38 
4] -31.74 42 
45 -30 46 
49 -31.99 50 
53 -44.79 54 
57 -34.38 58 





Ge] 
-30.69 
-30.5 
-30.22 
-30.01 
-30.12 
-30.91 
-33 
-38.25 
-56.52 
-35.68 
-30.99 
-30.09 
-33.39 
-66.67 
-32.29 


Degrees 
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G{9] 
-30,66 
-30.44 
-30.15 

-30 
-30.24 
-31.27 
-33.88 
-40.84 
-45.59 
-34.03 
-30.45 
-30.43 
-35.47 
-43.27 
-30.94 


Degrees 





G8] 
-30.62 
-30.36 
-30.09 
-30.01 

-30.4 
-31.73 
-34.98 
-44.86 
-40.87 
-32.74 
-30.12 
-31.03 
-38.72 

-37.6 
-30.19 








-30.01 
“44.5 
-21.67 
4732 
-6.62 
-3.04 
0.91 
-0.04 
-0.33 
-1.81 
-4.63 
OAL 
-16.21 
-30.67 
-31.96 
-30.94 
-43.27 
-35.47 
30.43 
-30.45 
-34.03 
-45.59 
-40.84 
-33.88 
31.27 
-30.24 
-30 
-30.15 
-30.44 
-30.66 
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-35.18 
-25.42 
-14.08 
-7.8 
-3.78 
-1.32 
-0.15 
-0.15 
-1.32 
-3.78 
-7.8 
-14.08 
-25.42 
-35.18 
-30.19 
-37.6 
-38.72 
-31.03 
-30.12 
-32.74 
-40.87 
-44.86 
-34.98 
-31.73 
-30.4 
-30.01 
-30.09 
-30.36 
-30.62 
-30.72 




















ARRAY CHARACTERISTICS 


Frew Hz): 2e00 
Celeritu €fesd: 42860 
* elements: is 
Spacinyw Cf td: e.36 
Wavetenevth (ft>: 2.43 
Look fingwle Cdew>: 98.0 
Window: TAYLOR 
DirectivitucdB>: 9,241 
2 AB hean width: 20.4 


FNSUREF PRINTFR FS RFADY 


PRESS ANY KEY TO PRINT 


Li 


Array Characteristics 


a 
a 
s 
~ 
a 

e 
- 

* 
= 
~ 
v 
9 





Frequency (Hz): 2000 
Sound Speed( fps): 4860 
# Elements: 15 
Element Spacing(ft): 0.6 
Wavelength (ft): 2.4 
Directivity (dB): 0.581 


Beamwidth (degrees): 20.44 


Element Amplitude Weight 
1 


2.457965 
4.369205 
6.542757 
8.704096 
10.54437 
11.78125 
6.10859 
11.78125 
10.54437 
8.704096 
6.542757 
4.369205 
2.457965 
1 
G[8] Degrees 
-46.22 2 
-46.67 6 
-47.94 10 
-50.79 14 
-58.75 18 
-57.95 22 
-47.36 26 
-42.51 30 
-40 34 
-39.6 38 
-42.64 42 
-70 46 
-39.15 50 
-33.81 54 
-34.64 58 
-41.26 62 


G{é] 

-46.27 
-46.9 

-48.46 
-51.98 
-64.23 
-54.1 

-45.86 
41.7 

-39.68 
-39.92 
-44.56 
51.44 
“47.2 

-33.37 
-36.64 
-33.45 


ARRAY BEAM PATTERN 


Degrees 
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Theta Cdewrees) 


G[9] 
-46.36 
47.18 
-49.08 
-53.54 

-70 
-51.32 
-44.58 
-41.01 

-39.5 
-40.49 
-47.61 
-45.39 
-35.71 
-33.29 
-40.82 
-28.72 


Degrees 





G6] 
-46.5 
-47.52 
-49.84 
-55.64 
-64.48 
-49.15 
-43.47 
-40.44 
-39.47 
-41.35 
-53.44 
-41.74 
-34.6 
-33.66 
-57.16 
-25.19 








-22.33 
-14.17 
-8.75 
-4.93 
-2.31 
-0.7 
-0.03 
-0.25 
-1.39 
-3.49 
-6.68 
-11.22 
-17.76 
-28.72 
-40.82 
-33,29 
-35.71 
-45.39 
-47.61 
-40.49 
-39.5 
-41.01 
-44.58 
-51,32 
-70 
-53.54 
-49.08 
-47.18 
-46.36 
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-17.76 
-11.22 
-6.68 
-3.49 
-1.39 
-0.25 
-0.03 
-0.7 
-2.31 
-4.93 
-8.75 
-14.17 
-22.33 
-41.26 
-34.64 
-33.81 
-39.15 
-70 
-42.64 
-39.6 
-40 
-42.51 
-47.36 
-57.95 
-58.75 
-50.79 
-47.94 
-46.67 
-46.22 





-15.87 
-9.93 
-5.77 
-2.87 
-1.02 
-0.11 
-0.11 
-1.02 
-2.87 
-5.77 
-9.93 

-15.87 

-25.19 

-57.16 

-33.66 
-34.6 

41.74 

-53.44 

41.35 

-39.47 

-40.44 

-43.47 

-49.15 

-64.48 

-55.64 

-49.84 

-47.52 
-46.5 
-46.2 











Array Characteristics 









ARRAY CHARACTERISTICS ARRAY BEAM PATTERN 


Freda (Wz): Jeoe 
Celerityu ¢€fps>: a360 
e elements: is 
spacing Cf t>: @.30 
Wavelength (ft): 8.97 
Look Angie Cdeg>: 939.9 
Window: USER 
Directivity(aB?: 13.8190 
2 AB heam width: 7.4 


ENSURE PRINTER IS READY 


G(Theta? (db) 


PRESS ANY KEY TO PRINT 


Theta (degrees) 





Frequency (Hz): $000 
Sound Speed( fps): 4860 
# Elements: ie) 
Element Spacing(ft): 0.6 
Wavelength (ft): 1 
Directivity (dB): 0.581 


Beamwidth (degrees): 7.12 


Element Amplitude Weight 


0.281302 
0.326266 
0.474929 
0.630153 
0.775678 
0.894747 
0.972813 
0.5 
0.972813 
0.894747 
0.775678 
0.630153 
0.474929 
0.326266 
0.281302 
G[9]} Degrees G[@] Degrees G[8] Degrees G8] 
-42.32 2 -41.9 3 -41.23 4 -40.39 
-39.42 6 -38.38 7 -37.31 8 -36.26 
-35.24 10 -34.27 Mi -33.37 12 -32.55 
-31.82 14 -31.2 15 -30.69 16 -30.31 
-30.08 18 -30.01 19 -30.13 20 -30.48 
-31.12 22 -32.14 23 -33.67 24 -36.02 
-39.92 26 -48.53 27 -50.46 28 -40.15 
-35.72 30 -33.06 31 -31.38 32 -30.39 
-30.01 34 -30.22 35 -31.12 36 -32.91 
-36.19 38 -43.18 39 -52.28 40 -38.58 
-33.9 42 -31.45 43 -30.25 44 -30.04 
-30.84 46 -32.91 47 -37.25 48 -50.5 
-41.61 50 -34,57 51 -31.48 52 -30.16 
-30.16 54 -34.57 55 -35.06 56 -44,36 
-43.49 58 -34.66 59 -31.28 60 -30.04 
li) 























Appendix C Digitized Hydrophone Output 


Pogea \14 Hour 


(32 Reproduced From 
Best Available Copy 
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3.1 Appendix Cl Unaveraged, Untapered, Hydrophone Data 











Hydrophone 1 with scattering plate Hydrophone | without scattering plate 











Hydrophone 2 with scattering plate Hydrophone 2 without scattering p!ate 





CROITIZED MYDROPHONE OUTPUT 








Hydrophone 3 with s -attering plate Hydrophone 3 without scattering plate 
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Hydrophone 4 with scattering plate 


DIGITIZED HYOROPHONE OUTPUT 


Hydrophone 5 with scattering plate 





CANITWZED HYOROPHONE OUTPUT 


Hydrophone 6 with scattering plate 








Hydrophone 4 without scattering plate 





OATA 








Hydrophone § without scattering plate 





ORR TUED MYOROPHONE OUTPUT 





Hydrophone 6 without scattering plate 
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Hydrophone 7 with scattering plate Hydrophone 7 without scattering plate 





CIGITIZED HYDROPHONE OUTPUT 
DATED HYDROPHONE OUTPUT 





Hydrophone 8 with scattering plate Hydrophone 8 without scattering plate 














Hydrophone 9 with scattering plate Hydrophone 9 without scattering plate 
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ORGITIZED HYOROPHONE OUTPUT 








Hydrophone 10 with scattering plate 











Hydrophone 11 with scattering plate 











Hydrophone 12 with scattering plate 
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Hydrophone 10 without scattering plate 





Hydrophone 11 without scattering plate 








Hydrophone 12 without scattering plate 























Hydrophone 13 with scattering plate 











Hydrophone 14 with scattering plate 











Hydrophone 15 with scattering plate 





Hydrophone 13 without scattering plate 





OMB TED HYOROPHONE OUTPUT 








Hydrophone 14 without scattering plate 











Hydrophone 15 without scattering plate 
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HYOROPHNE 
$id. sg ad 








Hydrophone 1 with scattering plate Hydrophone 1 without scattering plate 








Hydrophone 2 with scattering plate Hydrophone 2 without scattering plate 





OIAITZED NYDROPHONE OUTPUT 











Hydrophone 3 with scattering plate Hydrophone 3 without scattering plate 








Hydrophone 4 with scattering plate 


DIITZED HYDROPHONE OUTPUT 


Hydrophone 5 with scattering plate 





Hydrophone 6 with scattering plate 
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Hydrophone 4 without scattering plate 





OMB TED HYDROPMONE OUTPUT 
biissnees 





Hydrophone 5 without scattering plate 








ORR TIEED HYDROFHOME GUTRUT 
fide. aaay 





Hydrophone 6 without scattering plate 











Hydrophone 7 with scattering plate 


$e aa. 49 8 


DRT ZED » 





8 
RB 
g 


Hydrophone 8 with scattering plate 


; 
I: 


Hydrophone 9 with scattering plate 
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Hydrophone 7 without scattering plate 





Ox TEED HYDROPMONE CUTRUT 
Seis. saad 





Hydrophone 8 without scattering plate 


(ATED +vOROeoRE OUTPUT 
fidecsa 





Hydrophone 9 without scattering plate 

















QuTPYT 


CENITIZED HYOROPHONE OUTPUT 
fidecsaal 











Hydrophone 10 with scattering plate Hydrophone 10 without scattering plate 














E 1900 z 100 - 
1000 8 97000 + 
b _ 
; 1900 3 seco | 
pe 
° a Vtech al -m 
Hydrophone 11 with scattering plate Hydrophone 11 without scattering plate 





SSQITZED HYOROPHONE OUTPUT 
Gegec.a483 











iy 
g} 
g 


Hydrophone 12 with scattering plate Hydrophone 12 without scattering plate 
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Hydrophone 13 with scattering plate 








‘2000 
5 100 
f£ 
Bx 
wo 
° 
ro 
ts 
190 
Foe 
° 7 oa be 


Hydrophone 14 with scattering plate 


ROITIZED 1 VOROPHINE OUTPUT 








Hydrophone 15 with scattering plate 










OSE TUED HYOROMOME GCUTRUT 
Petes aead 


Hydrophone 13 without scattering plate 





biis.svtat 


Hydrophone 14 without scattering plate 









oureut 


YDADEHORE 
figs. adad 


Hydrophone 15 without scattering plate 


3.3 Appendix C3 Chebyshev Weighted Hydrophone Data (-30 dB) 


QATA 


oanzeD 


Hydrophone | with scattering plate 


DATA 


$988.99 9 3 





Hydrophone 2 with scattering plate 


QATA 
EY 





Once TED 





Hydrophone 3 with scattering plate 








bere had 





Hydrophone 1 without scattering plate 


Hydrophone 2 without scattering plate 





ONG TIZED MYDRDEHORE CUTRUT 
e 





Hydrophone 3 without scattering plate 
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OATA 





Hydrophone 4 with scattering plate 


QATA 


OaiTzED 


Hydrophone 5 with scattering plate 


DATA 





Hydrophone 6 with scattering plate 
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CAs TEED HYDROPHONE OUTRFUT 
Sigs cs gad 


Hydrophone 4 without scattering plate 





OKB TIED HYDROPMOS OUTPUT 





Hydrophone 5 without scattering plate 





ORE TEED HYORDFHOME OUTRUT 
hides. s dad 








Hydrophone 6 without scattering plate 




















DATA 


oan 
OXB TUE HYDROPHONE OUTPUT 





Hydrophone 7 with scattering plate Hydrophone 7 without scattering plate 








Hydrophone 8 with scattering plate Hydrophone 8 without scattering plate 





OATA 





OR nZED 


Ora TED HYDROPHORE OUTRFUT 
Phin. si ad 











Hydrophone 9 with scattering plate Hydrophone 9 without scattering plate 















ourpur 


iis. shit 


DATA 





QamnzeD 
omsrzen 








Hydrophone 10 with scattering plate Hydrophone 10 without scattering plate 





DRKTIZED HYDROPHONE DATA 





Hydrophone 11 with scattering plate Hydrophone 11 without scattering plate 








aara 





MYDROPHONME 
Pidt.rg at 


} 











Hydrophone 12 with scattering plate Hydrophone 12 without scattering plate 
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DAs TED MYOROPHONE DATA 








Hydrophone 13 with scattering plate 


ORAITIZED HYOROPHONE OUTPUT 


Hydrophone 14 with scattering plate 


DAGITHED HYOROPHONE OUTPUT 








Hydrophone 15 with scattering plate 
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Cad TEED HYOROFHONE GUTPUT 
fiis.ngad 








Hydrophone 13 without scattering plate 








Hydrophone 14 without scattering plate 








OATES +yOROPNOME OUTPUT 
iis .ntas 


Hydrophone 15 without scattering plate 








Appendix D Detecting the Scattered Signal 
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Hydrophone 1 


ie | 
me i 
| etme 








° xd ‘ain atic py 
Hydrophone 3 
ro 
os | 
Bs | 
| 
iz | 
' * TE (BECONDS) = es 
Hydrophone 5 
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Dear TEED HYOROPHONE OUTPUT 
figs. a8 83 





Hydrophone 2 
° a Said ccna ma a 
Hydrophone 4 





| 


OMITTED HYOROPHONE GU’ PUT 
Shes. a8 es 





Hydrophone 6 











L 





Gieecegag 


ALNO BeOrGdOMOAH CAT NO 


aon es ge 





Hydrophone 8 


Hydrophone 7 








gop eo ees g 


DLL BOP OMOQAH C1 WO 








TEL 


ANALOG INOMGOUGAH CGIZLUEND 





Hydrophone 10 


Hydrophone 9 











EPPe REE 


AD ELLAO DeOhSOWOAH GED 0 





Hydrophone 12 


Hydrophone 11 
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§$3° $3 8 
i ae 
. 


ANGLAIS Be POWOAH CEP INO 











Hydrophone 14 


Hydrophone 13 











Gite atiy- 


ADELNO JNOmeOeDAx CELLU 


Hydrophone 15 
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4.2 Appendix D2 Untapered Array Beamforming 


1.0x10* 


i 

z 0.5x10* 
= 

3 

Ww 

Zz 

fs) 

= 

a 

Ss) 0 
ac 

a 

> 

= 

Q 

8 -0.5x10* a 
9 

a 





-1.0x10* 
0 


TIME (SECONDS) 


Array response with plate in the water. 


1.0x10* 


0.5x10* 


-0.5x10* 


DIGITIZED HYDROPHONE OUTPUT 


-1.0x10* 





-1.5x10* 
0 


TIME (SECONDS) 


Array response with plate removed. 


133 














4.3 Appendix D3 Chebyshev Array Beamforming 


DIGITIZED HYDROPHONE OUTPUT 


DIGITIZED HYDROPHONE OUTPUT 


-2500 


-7500 


-4000 


7500 


2500 





TIME (SECONDS) 


Array response with plate in the water. 


8000 


4000 











02 06 


TIME (SECONDS) 


Array response with plate removed. 
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VOLTS?/Hz 


VOLTS*Hz 


.20 —- 





15 + 








0 2000 4000 6000 
Frequency (Hz) 


PSD omnidirectional hydrophone (#10) with plate in water. 


.20 





0 2000 4000 6000 
Frequency (Hz) 


PSD omnidirectional hydrophone (#10) with plate removed. 
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VOLTS*/Hz 


VOLTS?/Hz 











2000 4000 6000 
Frequency (Hz) 


PSD RECT window with plate in water. 


Frequency (Hz) 


PSD RECT window with plate removed. 
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1.25 a 





VOLTS*Hz 








0 2000 4000 6000 
Frequency (Hz) 


PSD Chebyshev Array with plate in water. 


1.6 


VOLTS*Hz 
° a @ x) 


PSD Chebyshev Array with plate removed. 
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